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_  '(U^This  report  presents  results  of  a  study  of  combustion  • D 

>rocesses  in  solid  propellant  cracks.  As  might  be  expoeted, 
mder  moderate  chamber  pressurization  rates  (<Hratm/s),  the 
:heoretical  predictions,  as  well  as  the  experimental  observations, 
Lndicate  that  the  ignition  front  propagates  from  the  entrance 
jfthecrac^  to  the  tip.  However,  under  rapid  chamber  pressurization 
rates  OKr atm/s  or  higher) ,  the  tip  region  of  the  crack  was 
>bserved  to  ignite  before  the  arrival  of  the  convective 
Ignition  front.  (Ignition  is  defined  here  as  the  onset  of 
•mission  of  luminous  light  from  the  propellant  surface  with 
some  material  loss.  A  theoretical  model  has  been  developed 
:o  explain  the  tip  ignition  phenomena.  The  model  considers: 
i  one-dimensional  transient  heat  conduction  equation  for  the 
solid  phase *  and  one-dimensional,  unsteady  mass  and  energy 
ronservation  equations  for  the  gas  phase  near  the  crack  tip. 

Joth  experimental  and  theoretical  results  indicate  that  the 
Ignition  delay  time  decreases  as  the  pressurization  rate  is 
Increased. yrfheoretically  calculated  ignition  delay  times  are 
in  good  agreement  with  the  experimental  data.  Based  upon 
:his  agreement,  ignition  near  the  crack-tip  region  is 
considered  to  be  caused  by  the  enhanced  turbulent  transport 
)f  energy  in  the  gas  phase,  which  is  driven  by  strong 
:ompression  waves. 

(U)  The  Crack  Combustion  Code  (CCC)  has  been  combined 
/ith  the  consolidated  Nonlinear  Finite-Element  Analysis 
Jrogram  (NFAP) .  Initial  runs  of  the  combined  program 
(CCC/NFAP)  indicate  that  the  structural  deformation  ind 
;ombustion  phenomena  are  interrelated.  The  combustion 
ihenomena  inside  the  crack  may  either  enhance  or  retard, 
lepending  upon  the  effect  of  mechanical  deformation  on  the 
5ap  width  along  the  crack.  The  program  has  the  option  of 
ising  either  the  linear  elastic  or  linear  viscoelastic  material 
nodel.  As  expected,  results  show  that  the  viscoelastic  model 
Ls  more  realistic  for  the  solid  propellant. 

(U)  Test  firings  were  conducted  to  determine  the  one- 
limensionelity  of  the  flame  front.  It  was  observed  that, 
sxcept  near  the  entrance  region,  the  one-dimensional 
issumpticn  for  the  flame  propagation  ls  quite  appropriate. 

Experiments  were  continued  to  determine  the  effect  of 
iropellant  crack  geometry  on  flame  propagation  inside  a 
:rack. 

(U)  A  complete  set  of  dimensionless  parameters  which 
govern  the  flame-spreading  and  combustion  phenomena  inside 
i  solid  propellant  crack  were  obtained. 
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INTRODUCTION 


A  very  important  characteristic  of  a  gas  permeable  solid  propellant 
is  that  it  undergoes  the  transition  from  normal  burning  to  detonation 
under  certain  operating  conditions.  A  transition  from  def lagration- 
to-detonation  (DDT)  may  not  only  cause  malfunction  or  failure  of  a  solid 
propellant  motor,  but  it  is  also  a  potential  source  of  hazard.  DDT 
usually  involves  very  high  flame  propagation  rates,  high  pressures,  and 
large  burning  surface  areas.  DDT  is  therefore  more  likely  to  occur  in 
gas  permeable  propellants,  such  as  granular  propellant  beds  or  propel¬ 
lant  grains  with  flaws  or  cracks  which  have  large  burning  surface  areas. 

Cracks  inside  solid  propellant  grains,  which  are  likely  to  be 
encountered  in  high-energy  propellants  with  high-solids  loading,  can 
allow  hot,  high  pressure  gases  to  penetrate  the  cavity.  The  combustion 
processes  inside  the  crack  can  produce  much  higher  pressure  than  the 
designed  maximum  operating  pressure.  If  the  local  pressure  rise  due  to 
the  gasification  is  sufficiently  rapid,  it  may  produce  strong  compres¬ 
sion  waves,  or  even  shock  waves  which  can  initiate  detonation.  Further¬ 
more,  for  DDT  in  damaged  propellant  grains  with  cracks,  crack  branching/ 
propagation  must  take  place  to  provide  the  additional  surface  area. 

Seme  experimental  and  theoretical  studies,  based  purely  on  solid  mech¬ 
anics  considerations,  have  been  conducted  to  predict  crack  propagation 
in  solid  propellant  grains.  However,  realistic  predictions  of  acceler¬ 
ative  convective  burning  and  crack  propagation  cannot  be  made  unless  a 
combustion  program  incorporating  a  structural  analysis  program  is  used. 

Convective  burning,  however,  is  not  the  only  mechanism  responsible 
for  producing  the  steep  pressure  gradients  that  lead  to  the  transition 
to  detonation.  There  are  other  mechanisms  which  may  also  be  important 
in  DDT  processes.  The  presence  of  numerous  ignited  spots  and  ablating 
surfaces  at  some  distance  from  the  combustion  zone  in  a  damaged  propel¬ 
lant  may  enhance  the  possibility  of  DDT.  For  example,  ignited  spots  can 
be  produced  by  reflection  of  shock-like  compression  waves  at  the  propel¬ 
lant  surface  located  in  the  region  ahead  of  the  combustion  zone.  A 
similar  phenomenon  was  observed  during  the  course  of  experimental  in¬ 
vestigations  of  ignition  and  f la, ue-sp reading  inside  a  solid  propellant 
crack.  As  noted  in  our  last  annral  report,*  under  rapid  pressurization 


1 Naval  Weapons  Center.  Convective  Burning  in  Isolated  Solid  Pro¬ 
pellant  Cracks ,  by  K.  K.  Kuo,  R.  L.  Kovalcin,  and  S.  J.  Ackman.  China 
Lake,  Calif.,  NWC,  February  1979.  (NWC  TP  6049,  publication  UNCLASSI¬ 
FIED.) 
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rates  (M.0  atm/s) ,  anamolous  ignition  was  observed  near  the  crack-tip 
region.  High-speed  (^15,000  pps)  morion  picture  films  useu  to  record 
the  experimental  observations  of  the  flame  propagation  inside  the  crack 
show  that,  at  high  chamber  pressurlzaton  rates,  the  tip  region  of  the 
crack  ignites  before  the  convective  ignition  front  propagates  from  the 
crack  entrance  to  the  tip.  Consequently,  under  these  operating  condi¬ 
tions,  two  flame  fronts  moving  in  opposite  directions  have  been  observed: 
one  propagates  from  the  crack  entrance  end  the  other  from  the  crack  tip. 

The  observed  ignition  phenomena  near  the  tip  region  could  be  one  of 
the  very  Important  mechanisms  responsible  for  the  generation  of  Ignited 
spots  in  a  region  ahead  of  the  combustion  zone.  Indeed,  these  ignited 
spots  can  contribute  significantly  to  the  DDT  process.  These  studies  in 
convective  burning  and  related  topics  are  essential  for  a  sound  physical 
understanding  of  the  phenomena  involved  in  the  DDT  process.  Moreover, 
solid  knowledge  of  the  ignition  phenomena  associated  with  strong  convec¬ 
tive  burning  is  an  important  step  towards  the  goal  of  modeling  a  more 
complex  DDT  phenomenon.  The  main  objectives  of  the  work  performed 
during  this  reporting  period  are  as  follows. 

1.  To  identify  the  possible  mechanism  for  ignition  near  the  crack 
tip  region,  and  to  develop  an  analytical  model  for  its  prediction. 

2.  To  consolidate  the  Nonlinear  Finite-Element  Analysis  Frogram 
(NFAP)  and  combine  it  with  the  Crack  Combustion  Code  (CCC)  . 

3.  To  experimentally  verify  the  one-dimensional  assumption  of  the 
flame  front. 

4.  To  determine  a  complete  set  of  dimensionless  parameters  that 
govern  the  combustion  phenomena  inside  a  solid  propellant  crack,  for  the 
combined  problem  of  aerothermochemistry  and  solid  mechanics. 

5.  To  continue  experimental  firings  of  single  pore  propellants 
with  narrow  gap-widths,  for  the  purpose  of  establishing  a  data  base  and 
for  further  verification  of  the  theoretical  model. 


BACKGROUND 

A  considerable  amount  of  work  in  the  area  of  solid  propellant 
crack  combustion  has  been  conducted  at  the  Pennsylvania  State  Univer¬ 
sity.1-3  Substantial  progress  has  been  made  in  both  experimental  and 

2K.  K.  Kuo,  A.  T.  Chen,  and  T.  R.  Davis.  "Convective  Burning  in 
Solid-Propellant  Cracks,"  AIAA  Journal,  Vol.  16,  No.  6  (June  1978),  pp. 
600-607 . 

3 Naval  Weapons  Center.  Transient  Combustion  in  Soli/1  Propellant 
Cracks ,  by  K.  K.  Kuo,  D.  R.  McClure,  A.  T.  Chun,  and  F.  C.  Lucas.  China 
Lake,  Calif.,  NWC,  October  1977.  (NWC  TP  5943,  publication  UNCLASSIFIED.) 
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theoretical  areas  towards  a  better  understanding  of  the  convective 
burning  inside  a  solid  propellant  crack.  Some  of  the  important  achieve¬ 
ments  are  as  follows: 

1.  A  quasi-one-dimenaional  theoretical  model  has  been  developed 
to  describe  the  transient  combustion  processes  in  a  propellant  crack. 

The  theoretical  model  can  be  used  to  predict  the  flame  propagation, 
pressurization,  and  gas  dynamic  phenomena  at  various  locations  along  the 
propellant  cavity. 

2.  A  computer  program,  based  upon  the  above  theoretical  model, 
has  been  developed  and  successfully  implemented.  This  program  can  be 
used  to  predict  the  flame  propagation  and  combustion  phenomena  inside 
the  propellant  crack. 

3.  A  laboratory-sized  combustion  chamber,  utilizing  current 
techniques  in  high-speed  photography  and  pressure  recording,  has  been 
developed  for  solid  propellant  crack  combustion  experiments  and  mode] 
validation. 

4.  In  order  to  simulate  the  conventional  pressurization  rates 
during  rocket  motor  ignition  transients,  a  reproducible  solid  propellant 
ignition  system  has  been  designed  and  fabricated. 

5.  Experimental  firings  with  single-pore  propellant  grain  have 
been  conducted  to  study  the  effects  of  crack  geometry,  chamber  pres¬ 
surization  rates,  propellant  type,  etc.,  on  flame  spreading  and  pressure 
variation  inside  the  cavity. 

6.  The  combustion  program  has  been  modified  to  consider  different 
ignition  criteria,  and  to  obtain  the  indepth  temperature  profile  of  the 
solid  propellant,  in  order  that  the  dynamic  burning  behavior  can  be 
studied. 


7.  An  extensive  literature  survey  on  the  combustion  phenomena 
inside  a  solid  propellant  crack,  and  its  relevance  to  the  process  ce 
DDT,  has  been  conducted. 


THEORETICAL  WORK 

MODELING  OF  THE  CRACK  TIP  IGNITION  PHENOMENA 
Introduction 


As  mentioned  earlier,  under  rapid  chamber  pressurization  conditions 
ignition  may  be  achieved  near  the  crack  tip  region — far  ahead  of  the 
convective  combustion  zone  in  an  isolated  crack.  It  is  conceivable  that 
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in  the  combustion  of  a  prefractured  or  granulated  propellant,  numerous 
Ignited  regions  can  be  generated  in  a  similar  manner  by  rapid  rates  of 
gasification  and  pressurization  in  the  combustion  zone.  These  ignited 
regions  can  propagate  further  to  produce  even  higher  gasification  rates 
in  the  combustion  chamber;,  or  in  an  actual  rocket  motor.  Flame-front 
propagation  rates  can  be  Increased  not  only  by  the  presence  of  these 
ignited  regions,  but  also  by  the  steepening  of  compression  or  shock 
waves  near  the  flame  front;  eventually,  the  waves  become  strong  enough 
to  initiate  a  final  transition  to  detonation.  Therefore,  a  fundamental 
understanding  of  the  basic  mechanisms  that  contribute  to  the  generation 
of  these  ignited  regions  would  be  tn  important  step  in  the  convective 
burning  and  DDT  study. 

The  observed  ignition  of  the  tip  region  of  the  crack  before  the 
arrival  of  the  convective  ignition  front  may  have  been  caused  by  one  or 
more  of  the  following  mechanisms: 

1.  Turbulent  diffusion  of  the  hot  product  gases  to  the  crack  tip 
region,  and  the  stagnation  heating  of  the  propellant  by  the  hot  gases  at 
the  crack  tip. 

2.  Compression/shock  wave  heating  of  the  gas  adjacent  to  the 
propellant  surface  and  the  stagnation  heating  of  the  propellant  by  the 
hot  gases  at  the  crack  tip. 

3.  Local  heat  release,  due  to  the  chemical  reaction  of  the  com¬ 
bustible  species  transported  to  the  crack  tip  by  the  initial  pressure 
gradient . 


4,  Radiative  heating  of  the  crack  tip. 

5.  Vaporization  of  the  fuel  binder  at  low  temperatures,  followed 
by  a  gas  phase  reaction  between  the  fuel  and  the  oxidizer  gases  in  the 
air  near  the  crack  tip  region. 

Ignition  near  the  tip  region  can  be  achieved  if  some  of  the  com¬ 
bustible  species  evaporate  from  the  propellant  surface  near  the  entrance 
of  the  crack  and  if  these  species  are  transported  by  the  initial  pressure 
gradient  to  the  crack-tip  region,  where  they  react  ai\d  ignite.  Under 
these  circumstances,  it  is  possible  to  observe  the  ignition  of  the  tip 
region  prior  to  the  arrival  of  the  convective  ignition  front.  In  order 
to  investigate  this  possibility,  an  inert  crack  was  formed  with  pronel- 
lant  only  at  the  tip.  Tests  conducted  with  this  type  of  inert  crack 
also  showed  ignition  at  the  tip,  thus  eliminating  the  third  mechanism 
noted  above . 

Some  preliminary  calculations  were  conducted  to  evaluate  the  effect 
of  radiative  heating.  Calculations  showed  that  radiative  heating  of  the 
propellant  at  the  tip  by  the  hot  gases  in  the  combustion  chamber  is  far 
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too  low  to  cause  ignition.  One  of  the  main  reasons  for  such  a  small 
amount  of  radiative  heating  Is  that  the  crack  length  to  hydraulic  diam¬ 
eter  ratio  is  on  the  order  of  80,  resulting  in  a  very  small  shape  fac¬ 
tor.  Therefore,  the  fourth  mechanism  mentioned  above  was  eliminated. 

The  fifth  mechanism,  i.e.,  the  gas  phase  reaction  mechanism  is 
quite  probable  for  propellants  using  fuel  with  low  vaporisation  temper¬ 
atures.  However,  the  gas  phase  reaction  may  not  be  very  appreciable  for 
propellants  with  fuel  binders  having  high  vaporization  temperatures 
(I.e.,  of  the  same  order  as  the  vaporization  temperature  of  the  oxidizer 
crystals).  Propellant  A,  for  which  all  tip  ignition  experiments  were 
conducted,  is  an  example  of  the  latter  type;  in  other  words,  its  fuel- 
binder  and  oxidizer  have  almost  identical  vaporization  temperatures. 

For  the  present  study,  therefore,  the  gas  phase  ignition  mechanism  may 
not  be  significant,  and,  hence,  can  be  neglected.  It  should  be  noted, 
however,  that  for  other  types  of  propellants  this  mechanism  may  be 
important. 

This  leaves  us  with  only  two  possible  mechanisms  that  should  be 
considered.  The  ignition  could  be  caused  by  a  combination  of  these  two 
mechanisms.  An  extensive  literature  survey  was  conducted  for  these  two 
heating  mechanisms.  The  heat  transfer  mechanism  near  the  closed  end  of 
a  long  isolated  pore,  with  a  time-varying  interface  temperature,  has  not 
been  studied  in  the  past.  Some  of  the  published  results  which  are 
related  to  the  problem  being  considered  are  discussed  in  the  following 
paragraphs . 

McAlevy,  Cowan  and  Summerf ield4 ' 5  studied  the  mechanism  of  ignition 
of  composite  solid  propellants  by  hot  gases  in  a  shock  tube.  The  objec¬ 
tives  of  their  study  were  to  measure  the  ignition  delays  of  test  propel¬ 
lants  heated  to  a  shock;  to  develop  a  theoretical  model  based  upon  the 
gas  phase  ignition  mechanism  in  order  to  calculate  the  ignition  delay; 
and  to  study  tue  effect  of  oxygen  concentration  on  the  ignition  delay. 
For  the  experimental  measurements,  the  specimen  was  flush-mounted  in  the 
end  wall  of  a  shock  tube.  The  propellant  was  heated  by  high  tempera¬ 
ture,  high  pressure,  stagnant  gas  behind  the  reflected  shock  wave. 
Photocells  were  used  to  detect  ignition.  In  the  theoretical  model,  it 
was  assumed  that  at  the  instant  of  shock  reflection  a  stagnant  column  of 
hot  gases  comes  in  contact  with  the  propellant  surface.  The  temperature 
of  this  column  of  gas  was  obtained  from  the  reflected  shock  relation¬ 
ships.  The  interface  temperature,  which  remains  constant  with  time, 


4R.  F.  McAlevy,  III,  P.  L.  Cowan,  and  M.  Summerf ield.  "The  Mech¬ 
anism  of  Ignition  of  Composite  Solid  Propellants,"  Solid  ~ rope Z lan t  Poc¬ 
ket  Research.  Volume  1,  Progress  in  Astronautics  and  Rocketry,  an 
American  Rocket  Society  series.  New  York,  Academic  Press,  1960,  pp. 
673-692. 

SR.  F.  McAlevy,  III.  "The  Ignition  Mechanism  of  Composite  Solid 
Propellants,"  Ph.D.  Thesis,  Princeton  University.  June  1960. 
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was  obtained  by  considering  two  uniform  temperature,  semi-infinite 
bodies  brought  in  contact  at  the  instant  of  shock  reflection.  Pyrolized 
fuel  vapor  generated  at  the  propellant  surface  propagates  upstream  at  a 
rate  governed  by  the  laws  of  mass  diffusion.  This  fuel  vapor  reacts 
with  the  gaseous  oxygen  and  ignition  takes  place  as  r  result. 

The  general  conclusions  of  McAlevy,  Cowan,  and  Summerfield  were 
that  the  observed  ignition  delay  is  the  time  required  for  sufficient 
fuel  to  react  with  the  gaseous  oxygen,  and  that  the  ignition  delay  time 
was  found  to  increase  with  the  decrease  in  oxygen  concentration  and 
pressure.  Since  this  was  the  first  studj  on  the  gas  phase  ignition 
model,  the  analytic  model  has  several  limitations.  The  assumption  of 
constant  propellant  surface  temperature  condition  is  not  very  realistic. 
Calculated  ignition  delay  time  based  on  their  model  was  an  order  of 
magnitude  higher  than  the  experimentally  observed  ignition  delay. 
Contrary  to  their  theoretical  predictions,  they  observed  a  longer  (two 
times)  ignition  delay  for  pure  fuel  than  for  composite  propellant  with 
the  same  fuel. 

Hermance.  Shinnar,  and  Summerfield  improved  the  gas  phase  ignition 
model  by  including  oxidizer  species  equations.  Later,  they  proposed  a 
more  realistic  model  for  pure  fuel,  which  included  the  effect  of  heat 
feedback  to  the  propellant  surface.  Unfortunately,  they  did  not  make 
comparisons  with  the  experimental  data  and,  therefore,  no  conclusions 
about  the  quantitative  prediction  capabilities  can  be  drawn.  Although 
the  effect  of  oxygen  concentration  (as  observed  by  McAlevy  ot  al.) 
supports  the  concept  of  gas  phase  ignition,  there  remains  a  large  number 
of  uncertainties  and  unknowns  which  could  strongly  influence  the  quanti¬ 
tative  predictions  of  ignition  delay.  Some  of  these  are:  detailed 
processes  involved  in  fuel  and  oxidizer  pyrolysis;  gasification  tempera¬ 
tures  of  various  ingredients  in  the  propellant;  the  effect  of  exothermic 
oxidation  at  the  propellant  surface;  and  the  catalytic  effect  of  oxi¬ 
dizer  gases  on  the  vaporizing  propellant  surface.  Because  of  these 
reasons  and  the  fact  that  the  studies  were  conducted  for  shock  tube 
conditions  where  the  effect  of  pressurization  rates  cannot  be  observed, 
the  results  of  the  gas  phase  ignition  model  cannot  be  used  for  the 
present  study. 


6C.  E.  Hermance,  R,  Shinnar,  and  M.  Summerfield.  "Ignition  of 
Hot,  Stagnant  Gas  Containing  an  Oxidizer,"  AIAA  J. ,  Vol.  3,  No.  9  (Sep¬ 
tember  1965),  pp.  1584-1592. 

7 - .  Ignition  of  an  Evaporating  Fuel  in  a  Hot  Oxidizing  Gas, 

Including  the  Effect  of  Heat  Feedback,"  Astvonautica  Acta,  Vol.  12,  No. 
2  (March- April  1966),  pp.  95-112. 
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Fay  and  Kemp  predicted  the  heat  transfer  rates  to  the  end  wall  of  a 
shock  tube.0  In  their  analysis,  the  end  wall  of  a  shock  tube  was 
assumed  to  be  suddenly  put  in  contact  with  the  quiescent  hot  gas.  The 
pressure  of  the  gas  behind  the  reflected  shock  was  determined  from  shock 
relationships.  Continuity  and  energy  equations  in  the  gas  phase  were 
used.  The  wall  was  assumed  to  remain  at  its  initial  temperature.  By 
using  an  integral  method,  a  similarity  type  of  solution  was  obtained  for 
the  heat  transfer  rate  to  the  end  wall.  A  closed  form  solution  was 
obtained,  which  predicted  the  heat  flux  as  a  function  of  shock  velocity 
and  time.  For  the  following  reasons,  results  obtained  from  the  study 
cannot  be  used  to  predict  the  heat  transfer  rate  in  the  present  study: 
a  shock  tube  type  of  study  cannot  take  into  account  the  effect  of  pres¬ 
surization  rates;  the  end-wall  temperature  is  ass  med  to  be  at  its 
initial  constant  temperature,  and,  therefore,  is  not  realistic  for 
ignition  studies;  and,  the  rate  of  heat  transfer  to  the  end  wall  can 
only  decrease  with  time  (q  ^1/^t)  in  the  analysis  (this  is  true  for  the 
Fay-Kemp  study,  but  is  contrary  to  the  physical  situation  of  the  present 
investigation) . 

Levy  and  Potter  made  some  transient  measurement  in  a  rarefaction 
wave  tube.8 9  Their  objectives  were  to  build  and  calibrate  a  heat  trans¬ 
fer  instrument,  to  perform  transient  heat  flux  measurements  in  a  rarefac¬ 
tion  wave  tube,  and  to  evaluate  the  significant  parameters  which  influ¬ 
ence  the  instrument  performance  in  the  measurement  of  the  transient  heat 
flux.  Using  a  heat  transfer  measuring  instrument,  thay  experimentally 
measured  the  heat  transfer  from  the  end  wall  for  various  pressure  ratios. 
In  their  theoretical  approach,  they  used  energy  and  continuity  equations 
in  a  control  volume  very  close  to  the  end  wall.  Temperature  of  the  end 
wall  was  assumed  to  be  constant,  and  the  gas  temperature  variation  was 
assumed  to  be  a  second  order  polynomial  of  the  distance  from  the  end 
wall.  The  expansion  process  was  assumed  to  be  isentropic.  From  their 
analysis,  they  found  that  was  a  function  of  a  dimen¬ 

sionless  time  tcin^tia^/Ltu^_F .  This  analysis  has  very  limited  use  for 

the  present  study  because  of  the  severe  simplifying  assumptions  and  the 
physical  process  of  depressurization. 


8J.  A.  Fay  and  N.  H.  Kemp.  "Theory  of  Heat  Transfer  to  a  Shock- 
Tube  End-Wall  from  an  Ionized  Monatomic  Gas,"  J.  Fluid  Meah . ,  Vol.  21 
(1965),  pp.  659-672. 

9M.  J.  Levy  and  J.  H.  Potter.  "Some  Transient  Measurement  in  a 
Rarefaction  Wave  Tube,"  J.  of  Engineering  for  Industry,  Trans.  ASME, 
Series  B,  Vol.  86,  November  1964,  pp.  365-370. 
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Chao  gave  an  alternative  approach10  to  the  same  problem  solved  by 
Levy  and  Potter.  Using  an  analytic  solution,  Chao  solved  the  energy  and 
continuity  equations  for  an  ideal  gas.  The  end  wall  was  assumed  to 
remain  at  its  initial  temperature.  Gas  density  was  assumed  to  be  in¬ 
dependent  of  the  axial  location.  By  using  a  coordinate  transformation, 
the  energy  equation  for  the  gas  phase  was  converted  into  the  form  of  a  • 
heat  equation.  A  closed  form  solution  for  the  dimensionless  heat  flux 
was  obtained.  The  predicted  heat  transfer  rate  was  in  close  agreement 
with  Levy  and  Potter's  data.  Once  again,  the  assumption  of  constant  end 
wall  temperature  and  the  physical  process  of  depressurization  makes  this 
analysis  of  very  limited  use  to  the  present  study. 

Theoretical  Model 


From  the  above  literature  survey,  it  is  clear  that  none  of  the 
existing  theoretical  models  car.  be  employed  to  study  the  current  igni¬ 
tion  phenomenon  observed.  Therefore,  a  theoretical  model  was  developed 
to  predict  the  ignition  delay  and  its  dependence  on  the  pressurization 
rate.  The  theoretical  model  considers  an  inert  crack — a  rectangular 
shaped  inert  channel  with  propellant  only  at  the  tip.  When  the  chamber 
is  pressurized,  hot  ignition  gases  penetrate  the  cavity  and  are  reflected 
from  the  solid  propellant  at  the  closed  end.  Ignition  delay  is  defined 
as  the  time  lag  between  the  arrival  of  the  pressure  front  at  the  tip  and 
the  incipient  ignition. 

As  soon  as  the  hot  gases  reach  the  closed  end  of  the  crack,  heat 
starts  to  flow  into  the  propellant  surface  and,  consequently,  the  sur¬ 
face  temperature  rises.  Ignition  is  defined  as  the  attainment  of  a 
critical  temperature  at  the  propellant  surface.  It  is  assumed  that  the 
solid  or  gas  phase  reaction  takes  place  before  the  onset  of  ignition. 
Here,  the  tacit  assumption  is  made  that  there  is  little  time  between  the 
attainment  of  a  critical  surface  temperature  and  the  reaction  between 
the  fuel  and  the  oxidizer  species  to  cause  ignition.  For  propellant  A, 
the  ignition  of  which  is  studied  here,  the  surface  temperature  ignition 
criteria  may  be  quite  suitable  since  its  fuel  binders  and  the  oxidizer 
have  an  almost  identical  vaporization  temperature.  It  should  be  pointed 
out  that  if  the  fuel  vaporization  temperature  is  much  lower  than  that  of 
the  oxidizer  crystals,  the  gas  phase  reaction  between  vaporized  fuel  and 
oxidizer  in  the  ambient  gas  may  be  important,  and  that  this  ignition 
criteria  may  not  be  valid. 


10B.  T.  Chao.  "End-Wall  Heat  Transfer  in  a  Rarefaction  Wave  Tube," 
J.  Heat  Transfer,  August  1  65,  pp.  349-352. 


10 


NWC  TP  6278 


In  the  theoretical  model,  the  propellant  is  considered  to  be  a  one¬ 
dimensional,  semi-infinite  solid  which  is  brought  in  contact  with  hot 
gases  at  the  Instant  of  the  arrival  of  the  pressure  front  ut  the  tip. 

The  hot  gases  cause  heat  transfer  to  the  propellant  surface.  A  schematic 
diagram  of  the  physical  model  considered  is  given  in  Figure  1. 


FIGURE  1.  Schematic  Diagram  of  the  Physical  Model 
for  the  Tip  Ignition  Study. 


The  solid  phase  heat  conduction  equation  is: 
?t  a2T 


The  initial  and  boundary  conditions  are: 

*  ■  °=  y°-x)  y 

x  -  :  yt.->  -  Tpl 

3T 

x " 0:  kP^  (t’0+>  * k M  <c’0~>  -  -v*> 


(1) 


(2) 
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where  $^(t)  Is  the  rate  of  heat  transfer  from  the  hot  gases  to  the 
propellant  surface  and  Is  unknown.  In  order  to  determine  the  heat  flux 
at  the  propellant  surface,  the  gas  phase  equations  have  to  he  considered. 

Assuming  the  bulk  flow  to  be  predominantly  one-dimensional,  the 
transient  mass  and  energy  conservation  equations  are  considered  for  the 
gas  phase.  Experimentally  obtained  pressure-time  trace  near  the  crack 
tip  region  is  used,  instead  of  solving  the  pressure  at  the  crack  tip 
through  the  use  of  momentum  equation.  The  region  of  interest  is  that 
adjacent  to  the  solid  propellant  surface.  The  governing  equations  can 
be  written  as  follows: 


Continuity  equation: 


+  IfiH  .  0 

3t  3x 


(3) 


Energy  equation: 


Hph)  a(puh) 
3t  <*x 


small 


3P  >r 
"  at  “  " 


small 


(k  il)  + 

dx  1  3x' 


dissipation 


(4) 


*<‘’°  >  -  P™..ur,d 


(5) 


In  order  to  simplify  the  analysis,  the  following  basic  assumptions  are 
made: 

1.  u  -|^  <<  for  the  region  of  interest,  i.e.,  locally  uniform 
pressure  near  the  propellant  surface. 

2.  The  viscous  dissipation,  <J>,  is  small. 

3.  The  specific  heat,  C^,  is  constant. 

4.  The  hot  gases  behave  as  ideal  gas. 

Using  these  assumptions  and  then  taking  the  Favre  Average11  of  the  gas 
phase  conservation  Eq.  (3)  and  (4),  we  get  (Appendix  A  summarizes  the 
basic  definitions  and  operation  rules  for  Favre  Averaging) : 


ilA.  Favre.  "Equations  des  Gar  Turbulents  Compressibles,"  J.  de 
Mechanique ,  Vol.  4,  No.  3  (1965),  pp.  361-390. 
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U  +  f-  (pu)  “  0 

at  ax 

■  aF  _  j__  ,)i__  ihx  4  1_  (V  il!l  _  ph"u") 

p  3T  +  pU  3x  at  3x  *Pr  3x'  +  3x  8x 


(6) 

(7) 


Considering  k  ~  to  be  negligible  in  comparison  with  ph"u"  and  defining 

oX 


,rt— n-  _  pt  ah 
-  oh  u  =  t— 


Prt  3x 


(8) 


The  energy  equation  can  be  written  in  the  following  fora: 

ah  .  ah  ap  .  i_  (iu  i£  +  ill  \ 

p  3r  +  pU  3x  “  at  3x  \Pr  3x  Pr  3x/ 


(9) 


The  initial  and  boundary  conditions  are 
I.C.  t  -  0  T(0,x)  *  T± 

B.C.l  x  -  T(t,-°°)  » 

B.C.2  x  -  0  T(t,0)  -  Tp(t,0)  (10) 


Solution  Method 

Integrating  the  continuity  equation  (Eq.  (6))  and  using  the  boundary 
condition,  at 

x  =  0  u  ■  0 


we  get 


pu 


/  -(ap/at)dx 
*o 


(11) 


Replacing  the  pQ  term  in  the  gas  phase  energy  equation  (Eq.  (9))  by 
Eq.  (11),  we  get: 


3P  _  3_  /y_  3h  ^t_  3h\ 
at  "  3x  y Pr  3x  Prt  3x/ 


(12) 


Now,  in  order  to  obtain  the  heat  flux  to  the  propellant  surface  and  the 
propellant  surface  temperature,  Eq.  (1)  and  (12),  together  with  the 
boundary  conditions  (Eq.  (2)  and  (10))  must  be  solved  simultaneously. 

An  attempt  was  made  to  solve  these  equations  analytically  to  obtain 
closed  form  solutions  for  the  propellant  surface  temperature.  The 
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partial  differential  equations  can  be  transformed,  in  the  transformed 
coordinates,  into  a  diffusion  type  of  equation  (see  Appendix  B ) .  Some 
difficulties  were  encountered  in  obtaining  analytical  solution,  as 
discussed  in  Appendix  2.  In  order  to  obtain  theoretical  results  to 
compare  with  experimental  data,  this  set  of  equations  is  solved  numer¬ 
ically. 

The  region  near  the  propellant  surface  is  the  most  interesting;  it 
is  here  that  large  temperature  gradients  exist.  For  mathematical  tract- 
ability,  the  infinite  x-region  is  transformed  into  a  finite  region, 
using  transformation  given  below.  The  transformed  coordinate  system 
also  makes  possible  a  finer  grid-sice  near  the  propellant  surface  in  the 
real  spatial  coordinate,  while  using  a  uniform  grid  in  the  transformed 
coordinate. 


For  the  gas  phase  region: 
-»  _<  x  <  0 

2  -  -[1  -  1/(1  -  Ax)] 
which  gives 

x  «  0  z  *  0 

X  *  — 00  z  **  —1 


3_ 

3x 


Ad  +  *)2  fj 


a2(i  a  ,)4  i- 


3z 


+  2a2<*  +  z>3  h 


(13) 


After  making  these  substitutions  into  the  gas  phase  energy  equation  (Eq. 
(12))  and  the  initial  boundary  conditions  (Eq.  (10)),  and  some  simpli¬ 
fication  as  well,  we  get 


3T 

f*  3p 

dz 

(1  + 

+ 1 

IlL.  +  !s_\  2A 2( 

1  +  zV 

9T 

at 

L/o  * 
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I.C.  t  -  0  T(0,z)  -  Tt 
B.C.l  z  -  -1  T(t,-1)  -  T 


B.C.2  z  -  0  T(t,0)  -  T  (t,0) 

P 

The  transformation  for  the  solid  phase  region  0  jC  x  £  » 

1 


(1  - 


1  +  AjX 


I 


(15) 


(16) 


which  gives 
x  *  0 


z  *  0 
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3  a  /,  ,2  a 
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2  *  Ai  d  ~  ^  (1  -  z)  3z 
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After  making  these  substitutions  into  Eq.  (1)  and  (2),  and  after  some 
simplifications,  we  get 
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1  P 


(t,0+) 
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(18) 


Numerical  Scheme 

The  coupled  set  of  particl  differential  equations  (namely,  the 
energy  equation  for  the  fluid  and  the  transient  heat  conduction  equation 
for  the  propellant),  was  solved  using  a  finite  difference  numerical 
technique.  The  most  important  solution  parameter  was  the  Interfacial 
temperature  as  a  function  of  time.  The  heat  flux  to  the  propellant 
surface  was  then  easily  computed  from  the  calculated  temperature  distri¬ 
bution.  The  finite  difference  equation  was  written  for  the  equations  in 
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transformed  coordinates.  As  mentioned  in  the  previous  section,  this 
transformation  made  possible  the  use  of  uniform  grid  spacing  in  *- 
coordinate,  thereby  bypassing  the  instability  associated  with  nonuniform 
grid  spacing,  which  would  have  been  necessary  for  the  x-coordinate.  The 
compression  of  the  axis  near  the  interface  was  controlled  by  the  param¬ 
eters  A  and  A..  Larger  values  of  A  and  A^  result  in  finer  grid  spacing 
near  the  Interface. 


The  two  partial  differential  equations  (Eq.  (14)  and  (17)),  to  be 
solved  are  of  the  form 


H 

at 


B(T,z) 


II 

as 


+  C(T,z)  $-1  +  l)(T,z) 

3* 


(19) 


The  backward  difference  method  was  used  to  approximate  the  time  deriva¬ 
tive  and  central  difference  for  spatial  derivatives.  The  equations  were 
linearized  by  evaluating  the  coefficients  B,  C,  and  D  at  an  intermediate 
time  step.  When 

Ttk  -  T(kAt,  iAx) 

Eq.  (1()  is  approximated  by 
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The  resulting  set  of  simultaneous  algebraic  equations  was  solved  by 
using  a  standard  tridiagonal  matrix  inversion.  At  each  time  step,  the 
properties  at  k  +  1/2  were  evaluated,  using  linear  interpolation  between 
k  and  k  +  1  levels.  Iteration  was  performed  at  every  time  step,  until 
the  solution  converged. 

Discussion  of  Results 

The  input  variables  used  in  a  typical  run  for  the  tip  ignition 
program  are  given  in  Table  1.  The  computation  is  terminated  when  the 
propellant  surface  temperature  reaches  the  specified  ignition  temper¬ 
ature  of  850*K.  Since  no  experimental  or  theoretical  results  are  avail¬ 
able  in  the  literature  for  turbulent  transport  properties  under  highly 
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transient  flow  conditions,  the  initial  selection  of  k  ^/k  was  somewhat 
arbitrary.  However,  it  was  selected  because  of  its  c?ose  agreement  with 
one  set  of  experimental  data.  In  all  later  computations,  the  same  value 
of  kgf^/k  was  used. 


TABLE  1.  Input  Variables  for  the 
Tip  Ignition  Computer  Program. 


M 

w 

m 

26.10  gf/g-mole 

Y 

- 

1.21 

P„. 

1.71  g/cm3 

Pt 

k 

m 

0.80  x  10  3  cal/s-cm-°K 

Pr  0 

a 

m 

0.0018  cm  7 s 

pr 

Tf 

- 

3000°K 

T. 

m 

850°K 

lgn 

Ti 

m 

298°K 

T  . 

m 

298  °K 

Pi 

k 

m 

2.085  x  10  ^  cal/s-cm-°K 

keff/k 

m 

300 

pi 

m 

1033.23  *«/cm 

Figure  2  shows  the  ignition  delay  as  a  function  of  the  pressuriza¬ 
tion  rate  at  the  crack  tip.  This  figure  also  shows  the  comparison  of 
predicted  and  measured  ignition  delays.  The  experimental  setup  and  data 
reduction  procedure  are  discussed  in  the  section  on  Experimental  Work. 

It  can  be  seen  that  the  ignition  delay  time  decreases  as  the  pressuri¬ 
zation  rate  is  increased.  The  slope  of  the  curve  is  larger  at  the  lower 
pressurization  rates.  The  nature  of  the  curve  may  help  to  explain  why, 
at  higher  pressurization  rates,  the  tip  region  of  a  real  propellant 
crack  may  ignite  prior  to  the  arrival  of  the  convective  ignition  front. 

At  lower  pressurization  rates,  the  ignition  delay  time  by  this  mechniism 
is  much  longer  (>1  ms);  therefore,  the  convective  ignition  front  will 
reach  the  tip  before  the  tip  region  ignites,  due  to  stagnation  point 
heat  transfer  associated  with  the  compression  waves  and  turbulent  energy 
transport  in  the  gas  phase.  However,  as  the  pressurization  rate  is 
increased,  the  ignition  delay  time  drops  sharply,  and  under  certain 
operating  conditions,  the  ignition  delay  time  may  be  smaller  than  the 
time  for  the  ignition  front  to  propagate  to  the  tip.  It  should  be  noted 
that  the  numerical  values  of  the  ignition  delay  times  and  pressurization 
rates  reported  here  should  not  be  arbitrarily  used  for  any  crack  geometry. 
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dP/dt  X  10"4.  atm/s 


FIGURE  2.  Comparison  of  Predicted  and 
Measure  1  Ignition  Delays. 


The  results  for  the  tip  ignition  tests  are  for  a  given  geometry  and, 
more  importantly  for  an  inert  crack;  the  numerical  values  of  the  result 
should,  therefore,  be  used  with  caution.  The  theoretical  predictions 
are  also  dependent  upon  the  value  of  An  increase  in  this  ratio 

will  cause  the  ignition  delay  time  versus  pressurization  rate  carve  to 
shift  downward,  whereas  a  smaller  ’Talue  of  ^  ,^/k  will  cause  it  to  shift, 
upward.  The  validity  of  the  theoretical  model  and  the  proposed  physical 
mechanism  is  supported  by  the  good  agreement  between  theoretical  and 
experimental  results 

Figure  3  presents  a  plot  of  propellant  surface  temperature  and  heat 
flux  to  the  propellant  as  a  function  of  time  for  a  given  pressurization 
rate.  The  propellent  surface  temperature  increases  almost  linearly  with 
time,  although  it  shows  some  curvature  at  the  beginning.  A  higher  pres¬ 
surization  rate  would  cause  the  slope  of  the  curve  to  steepen,  since, 
for  a  higher  pressurization  rate,  the  ignition  delay  will  be  shorter. 

The  heat  flux  to  the  propellant  surface  also  increases  with  time:  as 
time  increases,  the  temperature  gradient  at  the  propellant-gas  interface 
steepens.  The  slope  of  heat  flux  versus  time  curve  is  higher  in  the 
beginning  and  decreases  with  time.  At  a  higher  pressurization  rate,  the 
curve  will  shift  upward. 
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FIGURE  3.  Calculated  Heat  Flux  and 
Propellant  Surface  Temperature  vs. 
Time. 
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The  calculated  temperature  distribution  at  different  times  for  a 
given  pressurization  rate  is  plotted  in  Figure  4.  The  temperature  is 
plotted  as  a  function  of  the  transformed  dimensionless  coordinate,  z,  so 
that  the  detailed  variations  near  the  crack  tip  region  can  be  illustrated. 
The  temperature  profile  in  the  solid  propeMant  shows  that  the  gradient 
(5T/3z)  _+  [which  is  proportional  to  (3T/'3x)x-q+]  steepens  as  time 

increasisV  The  gradients  near  the  interface  appear  to  be  steeper  on  the 
propellant  side  than  on  the  fluid,  because  the  effective  (turbulent 
plus  molecular)  thermal  transport  coefficient  for  the  gas  is  higher  chan 
that  of  the  propellant.  Although  the  high  temperature  region  in  the  gas 
phase  seems  to  be  located  at  a  considerable  distance  from  the  solid-gas 
interface,  the  actual  physical  distance  to  the  high  temperature  region 
is  only  about  0.5  cm. 


LINKAGE  OF  CRACK  COMBUSTION  CODE  AND  NONLINEAR 
FINITE  ELEMENT  ANALYSIS  PROGR/M 

introduction 

As  mentioned  earlier,  for  DDT  in  damaged  propellant  grains  with 
cracks,  crack  branching/propagation  is  essential  to  provide  the  addi¬ 
tional  surface  area  for  anomalous  burning.  A  crack  combustion  program 
which  can  also  predict  stress  concentrations  at  the  tip  is  needed  for 
realistic  theoretical  predictions  of  the  accelerative  burning  inside  the 
propellant  crack,  crack  propagation,  etc.  Since  the  propellant  is 
deformable,  burning  inside  a  solid  propellant  crack  is  basically  a 
coupled  solid  mechanics  and  combustion  phenomen  .  The  combustion 
phenomenon  inside  the  crack  is  strongly  influenced  by  the  crack  geometry. 
The  geometry,  or  gap  width,  can  be  altered  for  the  following  two  reasons: 

1.  Mass  loss  due  to  gasification  or  ignition  of  the  propellant 
surface  along  the  crack  during  the  combustion  process,  leading  to  a 
variation  in  the  geometry. 

2.  Mechanical  deformation  of  the  propellant  due  to:  (a)  the 
pressure  loading  on  the  iuterior  surface  of  the  crack;  (b)  the  pressure 
force  of  the  combustion  chamber  acting  on  the  propellant  grain  near  the 
crack  entrance;  and  (c)  stresses  in  the  propellant,  due  to  pressure  and 
temperature  gradients  in  the  grain,  can  also  cause  the  crack  geometry  to 
change  significantly. 

Both  the  burning  rate  and  the  mechanical  deformation  are  governed 
by  pressure  acting  on  the  crack  surfaces.  On  the  other  hand,  a  change 
in  gapwidth  or  geometry  will  cause  the  pressure  distribution  to  vary. 

This  pressure  distribution  along  the  crack  will  strongly  influence  the 
deformations  and  stress  concentration  at  the  crack  tip.  It  is  obvious, 
therefore,  that  pressure  and  crack  geometry  are  strongly  interdependent 
in  a  burning  solid  propellant  crack. 
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FIGURE  4.  Calculated  Temperature  Distribution 
at  Various  Times. 


In  order  to  predict  the  effect  of  mechanical  deformation  on  the 
ignition  phenomena  inside  the  crack,  and  at  the  same  time  the  prediction 
of  stress  concentration  at  the  crack  tip  due  to  a  realistic  pressure 
distribution  along  the  crack,  the  crack  combustion  program  was  merged 
with  a  well-established  structural  analysis  program.  Of  the  numerous 
finite-element  structural  analysis  program,  the  Nonlinear  Finite-Element 
Analysis  Program  (NFAP)  (a  modified  version  of  the  NONSAP)1-* was 
considered  most  suitable  for  this  study.  This  program  is  a  user-oriented 
program  which  can  be  easily  modified  according  to  need.  It  is  also 
capable  of  handling  geometrical  nonlinearities  caused  by  the  large 
deformations  which  are  typical  of  propellants  under  high  pressure. 

Another  important  feature  of  the  program  is  its  restart  capability.  The 
NFAP  was  modified  to  suit  the  needs  of  the  current  problem.  A  linear 
viscoelastic  material  model  was  added  to  the  program  in  order  to  treat 
the  mechanical  properties  of  the  propellants  more  realistically. 

1  2Un;.versity  of  Akron.  NFAP- — A  Nonlinear  Finite  Element  Analysis 

Program,  by  S.  Prachuktam  and  T.  Y.  Chang.  Department  of  Civil  Engineer¬ 
ing,  Akron,  Ohio,  October  1976.  (Report  No.  SE  76-3,  publication  UNCLAS¬ 
SIFIED.) 

13r.  J.  Bathe  and  E.  L.  Wilson.  "NONSAP — A  Nonlinear  Structural 
Analysis  Program,"  Nuoleav  Enainet rina  and  Design,  Vol.  29,  No.  2  (1974), 
pp.  266-293. 
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Consolidation  NFAP 


For  computational  storage  efficiency  of  the  combined  program,  the 
original  version  of  NFAP  was  modified  and  consolidated.  The  unnecessary 
portions  of  NFAP  were  deleted.  In  the  input  phase  of  NFAP,  the  sub¬ 
routine  for  bandwidth  reduction  was  deleted.  Bandwidth  reduction  is 
essential  for  a  structure  whose  geometry  is  very  intricate  ano  for  which 
it  is  not  manually  possible  to  number  the  nodes  in  order  to  obtain  the 
narrowest  bandwidth.  For  ouv  well-defined  crack  geometry,  it  was  felt 
that  the  bandwidth  reduction  subroutine  was  unnecessary  The  subroutine 
for  >lculating  three-dimensional  pressure  loading  was  axso  deleted. 

Since  the  propellant  geometry  being  considered  Is  two-dimensional,  the 
subroutines  for  handling  one-dimensional  and  three-dimensional  input 
element  data  have  also  been  deleted. 

In  the  assembly  phase  of  NFAP,  the  subroutines  for  the  assembly  of 
one-dimensional  and  three-dimensional  stiffness  matrices  have  also  been 
removed.  The  original  NFAP  had  several  two-dimensional  material  models 
which  were  considered  unnecessary,  and,  hence,  were  deleted.  Instead,  a 
two-dimensional  linear  viscoelastic  material  model  was  added.  This  will 
be  discussed  in  detail  in  the  next  section.  In  the  incremental  solution 
phase  of  NFAP,  once  again  the  one-dimensional  and  three-dimensional 
equilibrium  iteration  subroutines  have  been  deleted.  The  one-dimensional 
and  three-dimensional  stress  calculation  subroutines  have  also  been 
eliminated.  Deletion  of  these  subroutines  has  resulted  in  a  significant 
reduction  in  the  size  of  NFAP,  which,  in  turn,  results  both  in  a  reduc¬ 
tion  in  the  storage  space  and  in  computational  time  and  cost. 

Viscoelastic  Material  Model 


A  linear  viscoelastic  material  model  has  been  implemented  to  de¬ 
scribe  the  behavior  of  the  propellant.  The  material  behavior  of  the 
solid  propellant  is  assumed  to  be  viscoelastic  in  shear,  and  elastic  in 
bulk.  At  present,  only  the  isothermal  condition  is  considered,  and  the 
viscoelastic  response  is  assumed  to  be  thermorheologically  simple.  The 
stress  and  strain  tensors  are  defined  as  follows: 


where 


ij 


stress  deviators 
strain  deviators 
stress  components 
strain  components 


su  -  “y  - 1/3  hj  ”kk  <21> 

'ij  '  £ij  -  1/3  sij  €kk  <22) 

Then,  the  stress-strain  relations  with  zero  initial  conditions  are 
written  in  two  parts. 
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The  shear  behavior  is  governed  by 

s±.  =  ^  Gxa  -  S')  jp-  e±j  (t')  dt'  (23) 

where  G1  is  the  relaxation  modulus  in  shear,  and  C  is  a  "reduced  time" 
due  to  the  temperature  above  a  reference  state.  For  a  constant  temper¬ 
ature  field,  5  is  defined  by 

£  =  <KT)t  (24) 

where  <J>  is  a  shift  function  due  to  temperature  effect  on  the  viscoelas¬ 
tic  behavior  of  the  propellant.  For  most  viscoelastic  materials,  it  is 
usual  to  consider 

M 

C1  ■  8o  +  S  *«  e'Sn,t  <25) 

ra=l 


The  bulk  behavior  is  governed  by 


3kk  '  3K  Ekk 


where  K  is  the  bulk  modulus. 

For  time-dependent  structural  analysis,  it  i|+^|cessary  to  form  an 
incremental  stress-strain  relationship .  When  (  )  is  a  quantity  at 

time  t+At  and  (  )  at  time  t,  the  incremental  stresses,  Ao . . ,  and 
strains,  Ae„,  are  defined  by  ^ 

.  t+At  t 

A  a.  .  =  o .  .  -  a . . 

ij  ij 


A  t+At  t  (')1\ 

-  Eij  -  (27> 

Then,  from  Eq.  (21)  to  (27),  we  can  derive  the  incremental  stress-strain 
relations  in  matrix  form. 


where 


(Ao(  =  [DVEl{Ae[  -  \oq\ 

|Aa|  =  |Aa22,  Aa33,  At23#  Aon| 
)Ae|  =  !ae22,  Ae33,  Ay23,  Ae11| 


la  l  -  V  B  j  c!L,  C*  C*  | 

I  of  ^  m<m  22  m  33  m  23  m  Ilf 

m=i 
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[dve]  -  [dve] 


(K+2/3  A)  (K-1/3A)  0  (K-l/3  A) 

(K+2/3A)  0  (k-l/3  A) 


(K+2/3  A) 


Symmetric 


where 


A  =  gQ  +  ^  gm(l  -  e  8m^t)/(6!n4)At) 


B  =  1  -  e“3m*At  (32) 

m 

and  the  term  .  has  a  recursive  relationship,  i.e., 

m  ij 

, ,  -6  <t>At. 

rt  -3  4>At  _t-At  8ai'‘  e  m  ) 

mCij  ‘  e  "  „CiJ  + - Tm? -  eij  (j3> 

»  t  t-At 

eu  *  eu  -  cu  <34) 

From  the  virtual  work  principle  and  the  relationship  of  Eq.  (28), 
one  can  derive  the  finite-element  stiffness  equations  for  a  time  inter¬ 
val  [t,  t+ut]. 


[k]jAv)  -  {if|  - 


where  Av  represents  the  increment  of  nodal  displacement  vector,  and  Af 
the  increment  of  applied  nodal  force  vector.  The  element  stiffness 
matrix  is  given  by 


k  =  f[ B]T  [Dve]  lB]  dv 


and  fQ,  a  force  vector  due  to  the  viscoelastic  reponse  of  the  material 


fo  -/[bT]  W  dv 


where  B  is  a  strain-nodal  displacement  transformation  matrix. 


The  material  model  described  above  has  been  written  in  the  form  of 
computer  subroutines  and  these  subroutines  were  implemented  into  NFAF 
for  the  purpose  of  conducting  the  viscoelastic  analysis  of  solid  propel¬ 
lant.  A  total  of  three  new  subroutines  were  prepared: 
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ELT2D9:  Allocate  array  storages  for  calculating  stresses,  strains, 

and  material  matrix. 

IVEL29:  Initialize  the  arrays  relating  to  the  formation  of  the 

material  matrix. 

VSEL.29:  Calculate  incremental  stresses  and  strains  and  form  the 

viscoelastic  material  matrix. 

Fcr  the  validation  of  the  viscoelastic  material  model,  the  program  was 
tested  for  some  simple  cases  for  which  analytical  solutions  were  avail¬ 
able.  The  agreement  between  the  numerical  and  analytical  solutions  was 
excellent. 

For  actual  computations  of  the  combined  program,  the  material 
properties  were  obtained  from  Mr.  C.  F.  Price  of  NWC,  China  Lake,  Cal¬ 
ifornia.  They  ate  as  follows: 

K  *  700,000  psi 

Gx  -  212  +  654  exp(-0.095t)  psi  (38) 

Combination  of  NFAP  and  CCC 


Approach.  For  the  transient  burning  inside  a  solid  propellant 
crack,  many  interrelated  physical  phenomena  should  be  considered  in  the 
theoretical  modeling.  The  combination  of  convective  burning  and  struc¬ 
tural  mechanics  analyses  can  be  very  complex.  Some  of  the  interdepend¬ 
ent  processes  are: 

1.  Pressure  inside  the  crack  results  in  a  deformation  of  the 
propellant  crack  geometry. 

2.  During  combustion,  there  is  mass  loss  due  to  the  gasification 
of  the  propellant  surface  inside  the  crack,  and  this  changes  the  geom¬ 
etry  of  the  crack. 

3.  Deformation  of  the  propellant  crack  geometry  causes  the  pres¬ 
sure  to  change. 

4.  Pressure  wave  phenomena  is  interdependent  with  the  transient 
variation  in  the  crack  geometry. 

5.  Compression  of  the  propellant  may  alter  its  burning  charac¬ 
teristics. 

6.  The  temperature  gradients  involved  in  the  convective  burning 
may  cause  thermal  stresses  in  the  solid,  contributing  to  crack  propa¬ 
gation. 
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7.  The  structural  analysis  should  account  for  the  mass  loss  and 
the  effect  of  surface  ablation  on  the  stresses  in  the  solid. 

It  is  obvious  that  the  coupling  between  the  pressure  along  the  crack  and 
the  deformation  is  the  most  important,  and  that  it  could  influence  the 
combustion  phenomena  significantly.  Some  of  the  other  coupled  processes 
mentioned  above  are  of  secondary  importance  and  are  ignited  in  the 
present  analysis.  The  thermal  stresses  due  to  che  temperature  gradients 
are  considered  to  be  small  when  compared  with  the  mechanical  stresses 
caused  by  deformation  alone.  Variation  in  the  burning  rate  due  to 
compression  of  the  propellant  may  not  be  significant,  considering  the 
strong  dependence  of  the  burning  rate  on  pressure. 

Incorporation  of  the  mass  loss  of  the  propellant  into  the  struc¬ 
tural  analysis  may  be  important,  but  it  is  quite  involved,  requiring 
extensive  program  revisions.  However,  it  should  be  noted  that  for  the 
region  of  interest  of  our  analysis,  i.e.,  during  the  early  phase  of 
combustion  before  substantial  burning  takes  place,  ignoring  mass  loss  in 
the  structural  analysis  may  be  justified.  During  this  period,  the 
geometry  variation  due  to  surface  gasification  is  small  when  compared  to 
the  variation  due  to  mechanical  deformation,  which,  of  course,  is  ac¬ 
counted  for  in  the  program. 

In  the  present  solution  scheme,  the  pressure  and  the  deformation 
are  not  solved  simultaneously.  They  are  obtained  independently  at  a 
finite  time  interval  apart,  i.e.,  pressure  at  time  t^  is  used  as  an 
input  to  determine  the  deformation  at  time  t.+At.  In  other  words,  a 
time  lag  is  introduced  between  the  pressure  computations  obtained  from 
CCC,  and  the  deformation  computations  obtained  from  NFAP.  The  smaller 
the  time  lag,  the  more  accurate  the  solution  will  be.  A  schematic 
diagram  of  the  coupling  is  shown  in  Figure  5.  It  is  believed  that  the 
time  lag  between  the  pressure  and  deformation  computations  will  result 
in  only  negligible  errors.  In  fact,  it  may  even  be  justified  on  the 
grounds  that  the  material  would  require  some  finite  time  lag  to  respond 
to  an  applied  pressure  load.  Also,  these  time  lags  are  usually  smaller 
than  a  few  hundreds  of  the  total  event  time,  and,  therefore,  may  not 
have  a  very  significant  effect  on  the  overall  solution. 

Implementation  of  the  Combined  Program.  The  initial  implementation 
of  the  combined  program  included  such  considerations  as  input  and  output 
data,  an  overlay  structure  to  conserve  memory  core,  and  organization  of 
COMMON  storage  areas  to  transfer  necessary  data  between  the  two  pro¬ 
grams.  A  system  of  interfacing-subroutines  were  used  to  transfer  pro¬ 
gram  control  between  the  NFAP  routines  and  the  CCC  routines.  When 
coupling  the  two  programs,  the  approach  was  to  keep  them  as  separate  as 
possible,  for  the  following  two  reasons.  First,  if  modifications  to  CCC 
are  made,  program  modification  and  testing  can  be  made  for  CCC  alone, 
and  NFAP  and  CCC  can  be  recoupled  without  further  revisions.  This 
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ccc 

NFAP 

CCC/NFAP 

COMBINATION 


-  PRESSURE  DEPENDS  ON  TIME  HISTORY  AND 
DEFORMATION  OF  GAP  WIDTH 

-  MECHANICAL  DEFORMATION  IS  A  FUNCTION 
OF  PRESSURE  AND  TIME 

-  COUPLED  COMBUSTION  AND  SOLID 
MECHANICS  SOLUTION  FOR  ANY  SPECIFIED 
INCREMENTAL  TIME  STEP 


FIGURE  5.  Schematic  Diagram  of  CCC  and  NFAP  Linkage. 


improves  the  efficiency  of  any  program  testing,  and  allows  the  greatest 
flexibility  in  program  verification.  Similarly,  NFAP  can  be  tested 
alone,  with  a  coupling  to  CCC,  or  with  an  artificial  pressure  loading 
routine,  in  order  to  perform  solid  mechanics  studies.  A  second  reason 
for  maintaining  the  two  programs  separately  is  to  keep  NFAP  and  CCC  as 
close  to  their  original  forms  as  possible.  In  this  way,  others  can 
utilize  the  programming  procedures  developed,  substituting  their  own 
procedure  for  convective  burning  or  structural  analysis,  if  desired, 
without  formulating  a  new  linkage.  These  considerations  increase  the 
potential  usefulness  of  the  combined  program  by  keeping  it  general  and 
versatile . 

The  disadvantages  of  this  approach  are  minimal.  The  greatest 
disadvantage  of  keeping  NFAP  and  CCC  as  complete  as  possible  is  the  need 
for  increased  program  memory  size.  Both  are  large  programs,  with  large 
core  storage  requirements.  This  limitation,  however,  can  be  relaxed  by 
utilizing  a  memory  overlay  structure  to  conserve  memory  core  for  program 
use.  Since  the  two  halves  of  the  combined  NFAP/CCC  program  are  executed 
alternately,  only  one  half  needs  to  be  in  memory  at  any  given  instant 
during  the  execution  of  the  program.  That  is,  while  the  structural 
analysis  is  taking  place  during  the  execution  of  the  combined  program, 
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the  convective  burning  routine  is  not  immediately  needed,  so  the  space 
in  memory  is  overlaid  by  the  NFAP  routines.  Similarly,  when  the  con¬ 
vective  burning  analysis  is  taking  place,  the  NFAP  routines  not  in 
execution  at  that  time  are  replaced  in  core  by  the  CCC  routines.  The 
penalty  of  increased  execution  time  due  to  the  transfers  between  memory 
core  and  an  intermediate  input/output  device  (direct  access  magnetic 
disk)  is  small  when  compared  to  overall  execution  time.  The  overlay 
procedure  is  automatic  and  independent  of  the  program.  The  operating 
system  for  the  IBM  Model  370  computer  used  for  all  the  computations 
includes  an  overlay  feature  which  requires  only  an  outline  of  the  over¬ 
lay  structure  and  some  control  cards.  The  use  of  the  overlay  feature 
was  not  critical  for  any  of  the  analyses  performed  on  the  IBM  370  com¬ 
puter  at  Pennsylvania  State  University  because  of  the  large  memory 
capacity  (1000  k  bytes  of  storage  to  a  single  program)  available.  The 
core  size  specification  of  the  operating  system  allowed  maximum  avail¬ 
able  storage  and,  consequently,  the  overlay  feature  was  not  needed  for 
the  computations.  At  other  computer  installations,  however,  core  size 
limitations  are  usually  very  critical,  and  an  overlay  structure  will  be 
necessary. 

The  theoretical  coupling  considered  between  the  convective  burning 
and  structural  analysis  is  simplified,  as  previously  discussed,  to  the 
transfer  of  pressure  loading  and  crack  deformation.  Pressure  is  calcu¬ 
lated  by  the  CCC  portion  of  the  combined  program  at  each  nodal  point 
location  on  a  one-dimensional  grid  along  the  length  of  the  crack.  The 
convective  burning  analysis  of  the  crack  combustion  incorporates  the 
crack  geometry  variation  caused  by  both  mechanical  deformation  and  mass 
loss  through  gasification  of  the  propellant  surface.  Once  the  gas 
phase  equations  are  solved  and  pressures  along  the  crack  are  calculated 
for  a  particular  time  t,  program  control  is  transferred  to  the  NFAP 
portion  of  the  combined  program.  The  surface  elements  in  the  finite 
element  mesh  are  then  loaded  with  the  pressure  calculated  in  CCC.  NFAP 
updates  the  stiffness  matrix,  performs  a  static  analysis  at  time  t,  and 
calculates  the  propellant  deformation.  The  change  in  the  crack  width  at 
each  finite-difference  node  is  calculated  and  added  to  the  existing 
crack  width.  Control  is  transferred  back  to  CCC,  the  next  time  itera¬ 
tion  (t  +  At)  is  made,  and  control  is  again  transferred  to  NFAP.  In 
this  manner,  pressure  and  changes  in  crack  width  are  coupled.  The 
general  layout  of  the  linkage  is  shown  in  Figure  6. 

Some  modifications  to  the  procedure  described  above  are  necessary 
due  to  considerations  of  execution  time  efficiency.  The  time  required 
for  NFAP  to  complete  a  static  analysis  is  approximatly  an  order  of 
magnitude  higher  than  the  execution  time  for  one  CCC  time  iteration. 
Although  CCC  uses  an  implicit  solution  method,  the  time  must  be  small 
enough  to  give  accurate  results.  Performing  structure  analysis  for 
every  iteration  in  CCC  is  prohibitive  in  view  of  the  extremely  long 
computer  time  required.  An  option  of  the  combined  program  coupling  is 
to  specify  the  number  of  time  iterations  CCC  will  complete  before 
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FIGURE  6 


General  Layout  of  the  CCC  and  LFAP  Linkage 


transfer  is  made  to  NFAP  for  an  updating  of  the  crack  width.  When 
CCC  takes  more  than  one  time  step  between  crack-width  updating,  the 
change  in  crack  width  is  averaged  over  the  time  increments  before  addi¬ 
tion  to  the  existing  width.  For  example,  if  CCC  completes  ten  itera¬ 
tions  before  transferring  control  to  NFAP  for  an  updating  of  crack 
geometry,  one-tenth  of  the  change  in  deformation  returned  by  NFAP  will 
be  added  during  the  next  ten  iterations  in  CCC.  This  averaging  approx¬ 
imation  is  designed  to  smooth  out  any  step  changes  in  the  geometry  that 
would  affect  the  CCC  solution  convergence.  The  allowable  number  of 
x.. orations  in  CCC  between  transfers  to  NFAP  is,  in  general,  dependent  on 
the  time  increment,  chamber  pressurization  rate,  and  crack  width. 

The  general  stress-strain  equations  are  solved  using  a  plane-strain 
analysis.  This  is  compatible  with  the  experimental  test  configuration, 
in  which  the  propellant  is  confined  from  all  sides  except  the  entrance. 
The  tra^Tient  nature  of  the  pressure  loading  is  considered  by  using  a 
statif  ..aalysis  at  incremental  time  steps.  A  static  analysis  was  chosen 
because  of  the  prohibitive  computational  times  involved  and  the  complex 
natu’-'  of  dynamic  analysis.  However,  the  options  of  using  the  dynamic 
an  1  ■  Is  is  available. 

emulation  of  the  Finite  Element  Mesh.  The  analysis  considers  a 
two-d..^ensional  configuration,  which  is  compatible  with  the  crack  geom¬ 
etry  considered  in  the  convective  burning  model  in  the  crack  combustion 
coae.  TV-' -dimensional  continuum  elements  were  used  to  model  a  sample 
crack.  ,  cause  of  symmetry  about  the  crack  centerline,  only  one  half 
of  the  solid-propellant  crack  model  needs  to  be  considered.  Figure  7  is 
a  diagi  of  a  finite  element  grid  for  a  sample  problem.  The  element 
mesh  near  the  crack  entrance  was  slightly  modified.  The  square  corner 
at  the  entrance  of  the  crack  was  rounded  a  bit,  and  the  grid  was  refined 
somewhat  for  better  computational  accuracy.  Each  isoparametric,  quad¬ 
rilateral  element  contains  eight  nodes,  four  nodes  at  the  corners  and 
one  node  at  the  midpoint  of  each  side.  This  feature  of  NFAP  makes  it 
possible  to  reduce  the  number  of  elements  needed  to  produce  accurate 
results.  Equal  accuracy  may  thus  be  achieved  with  less  computational 
effort  than  is  necessary  with  four-node  elements.  During  the  initial 
testing  of  the  program,  a  simplified  element  grid  was  used,  with  only 
one  layer  of  elements  instead  of  the  present  configuration  of  four  layers 
in  the  axial  direction  of  the  crack.  Multiple  layers  are  used  in  order 
to  increase  the  accuracy  of  the  solution. 

The  finite-difference  nodes  used  in  the  CCC  one-dimensional  grid 
are  in  a  one-to-one  correspondence  with  the  finite-element  nodes  which 
form  the  propellant  surface  in  the  NFAP  two-dimensional  mesh.  This 
allows  one  to  transfer  the  deformations  obtained  from  NFAP  directly  to 
CCC.  The  model  grid  is  stressed  by  applying  pressure  loading  to  the 
appropriate  surfaces,  and  fixing  boundary  conditions  to  the  elements. 

The  transient  pressure  loading  for  the  individual  elements  along  the 
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FIGURE  7.  Finite  Element  Grid. 


crack  surface  is  determined  by  the  crack  combustion  routine.  The  load¬ 
ing  at  the  entrance  of  the  crack  is  determined  by  the  chamber  pressure, 
given  as  a  function  of  time.  The  schematic  diagram  of  the  pressure 
loading  and  boundary  conditions  of  the  model  grid  is  shown  in  Figure  8. 
All  surfaces  of  the  propellant  mesh  are  either  pressure  loaded  or  have 
fixed  boundaries.  The  propellant  at  the  tip  of  the  crack  is  fixed  in 
one  direction  only,  due  to  the  symmetry  of  the  mesh.  Along  the  surfaces 
of  the  propellant  in  contact  with  the  brass  mold,  the  boundaries  are 
assumed  fixed. 

Pressure  data  from  the  convective  burning  analysis  are  available  at 
the  finite-difference  nodes.  Pressure  input  to  NFAP,  however,  must  be 
in  the  form  of  uniform  pressure  loading  on  individual  element  surfaces. 
Every  element  surface  is  defined  by  three  nodes:  two  nodes  at  the 
endpoints  of  the  elements,  and  one  node  at  the  midpoint  of  the  surface. 
The  uniform  surface  loading  of  an  element  is  obtained  from  the  spatially- 
weighted  average  of  the  pressures  at  the  three  surface  nodes.  Although 
there  is  a  one-to-one  correspondence  between  the  finite-element  and 
finite-difference  nodes  along  the  surface  of  the  crack,  the  pressure 
profile  considered  by  NFAP  contains  discontinuities  at  each  element 
boundary  along  the  surface,  because  of  the  pressure  averaging  procedure. 

Pressure  loading  of  the  two-dimensional  elements  at  the  entrance 
and  the  tip  of  the  crack  also  involved  side  loading,  i.e.,  loading  in 
the  axial  direction  of  the  crack.  The  side  loading  of  the  surface 
elements  at  the  opening  of  the  crack  was  determined  by  the  chamber 
pressure  at  any  given  time.  Side  loading  at  the  tip  of  the  crack  was 
considered  the  same  as  pressure  loading  on  the  last  element  along  the 
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surface  of  the  crack.  Appendix  C  lists  the  input  format  and  Anpendix  D 
the  sample  data  for  the  combined  program. 


FIGURE  8.  Boundary  Condition  for  a  Test  Specimen. 


Discussion  of  Results.  Some  preliminary  computations  were  made 
with  the  combined  structural  analysis  and  convective  burning  program. 

The  results  obtained  from  the  combined  program  were  compared  with  those 
obtained  from  convective  burning  analysis  alone.  From  these  compari¬ 
sons,  it  was  possible  to  draw  some  general  conclusions.  A  more  detailed 
analysis  of  the  results  from  the  combined  program  is  being  performed 
under  the  second  phase  of  the  contract. 

The  crack  geometry  used  for  the  combined  NFAP/CCC  program  is  com¬ 
patible  with  the  dimensions  of  the  propellant  test  specimen  used  for 
experimental  test  firings.  The  primary  inputs  to  the  program  are  the 
chamber  pressurization  rate  and  the  crack  geometry.  For  preliminary 
test  runs,  the  initial  geometry  of  the  crack  has  a  uniform  rectangular- 
cross-scctional  area.  The  chamber  pressurization  rates  are  assumed  to 
be  constant,  i.e.,  the  chamber  pressure  is  a  ramp  function  of  time.  For 
a  very  low  chamber  pressurization  rate,  the  effect  of  the  mechanical 
deformation  is  relatively  insignificant  on  the  combustion  process  inside 
the  crack. 

The  chamber  pressurization  rate  is  5  x  10^  atm/s  for  the  result 
reported  here.  It  should  be  pointed  out  that  the  constant  pressuriza¬ 
tion  rate  used  is  quite  close  to  that  observed  experimentally,  especial¬ 
ly  uver  the  period,  of  interest.  In  an  actual  test  firing,  the  chamber 
pressure  rises  at  a  nearly  constant  rate,  remains  steady  for  a  long 
period  of  time,  and  then  falls.  The  calculation  is  continued  until  the 
pressure  reaches  the  peak  chamber  pressure. 
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Calculated  pressure  distributions  at  various  times,  from  the  con¬ 
vective  burning  analysis  alone,  are  shown  in  Figure  9.  As  time  pro¬ 
gresses,  the  chamber  pressure  increases,  hot  gases  enter  the  crack,  and 
pressure  waves  travel  along  the  length  of  the  crack.  As  the  hot  gases 
flow  over  the  propellant  surface,  the  surface  temperature  increases 
until  surface  ignition  occurs.  Propellant  mass  loss  through  gasifica¬ 
tion  causes  increased  pressure,  which  enhances  the  burning  rate.  The 
combustion  gases,  driven  by  the  pressure  gradient,  penetrate  deeper  into 
the  crack,  causing  the  ignition  front  to  propagate.  The  locus  of  the 
ignition  is  shown  by  the  dashed  line  on  the  same  figure.  At  232.5  us 
the  ignition  front  has  propagated  through  60%  of  the  crack  length. 
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FIGURE  9.  Calculated  Pressure  Distributions 
for  Various  Times  From  the  Crack  Combustion 
Code . 


Using  the  same  initial  geometry  and  pressurization  rate,  the  test 
case  was  run  again,  but  with  the  coupled  convective  burning  and  struc¬ 
tural  analysis  program  incorporating  the  linear  viscoelastic  model.  The 
pressure  distributions  and  the  locus  of  the  ignition  front  calculated 
from  the  combined  NFAP/CCC  program  is  shown  in  Figure  10.  The  general 
trends  of  pressure  with  respect  to  time  are  similar  to  those  from  the 
convective  burning  analysis  alone,  but  when  the  pressure  inside  the 
crack  increases  above  25  atm,  the  pressure  distributions  from  the  two 
analyses  begin  to  differ.  The  pressure  obtained  from  the  combined 
analysis  generally  increases  faster  than  the  pressure  calculated  by  the 
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FIGURE  10.  Calculated  Pressure  Distributions  for  Various  Times 
From  the  Combined  Crack  Combustion  and  Nonlinear  Finite-Element 
Analysis  Program. 


convective  burning  orogram  alone.  This  can  be  explained  by  an  examina¬ 
tion  of  the  port  area  variation  along  the  crack.  The  increased  pressure 
corresponds  to  a  reduction  in  the  area  of  the  flow  channel.  The  physi¬ 
cal  explanation  of  the  decrease  in  crack  width  in  the  combined  analysis 
is  based  on  the  solid  mechanics  consideration.  Because  of  the  pressure 
loading  on  the  elements  at  the  entrance  of  the  crack  and  the  constraints 
along  the  boundaries,  the  deformation  of  the  propellant  was  fairly 
complicated.  The  elements  at  the  entrance,  compressed  by  the  chamber 
pressure,  deform  inward,  pushing  the  propellant  into  the  crack.  Since 
the  chamber  pressure  increases  faster  than  the  pressure  inside  the  crack 
during  the  early  stage,  the  propellant  is  pushed  towards  the  lower 
pressure  region  inside  the  crack.  The  mechanical  deformation  of  the 
propellant  causes  the  narrowing  of  the  crack  width. 

As  time  progresses,  the  pressure  distribution  becomes  wavy  in 
shape,  causing  the  propellant  surface  also  to  become  wavy.  The  local¬ 
ized  pressure  rises  correspond  to  the  deformation  of  the  pronellant 
along  the  crack.  The  pressure  and  port  area  display  a  very  strong 
inverse  relation.  Propagation  of  the  localized  pressure  phenomena  along 
the  crack  toward  the  tip  can  be  seen  in  Figure  10.  At  232.5  us,  the 
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localised  pressure  phenomena,  evident  around  nodes  12  and  20,  are  in¬ 
creasing  in  magnitude.  The  locus  of  ignition  front  exhibits  a  slightly 
faster  rate  of  propagation  into  the  crack,  as  compared  to  that  shown  in 
Figure  9.  It  Is  interesting  to  note  that  the  result  obtained  from  the 
coupled  program  does  not  necessarily  give  lower  pressure  inside  the 
crack  or  a  lower  flame  propagation  rate. 


DETERMINATION  OF  COMPLETE  SET  OF 
DIMENSIONLESS  PARAMETERS 


The  nondime ns Iona Used  governing  equations  for  the  gas-dynamic  and 
heat-transfer  equations  were  reported  in  the  last  annual  report,2  and 
will  not  be  repeated  here.  The  nondimens ional  equations  for  the  solid 
mechanics  portion,  and  the  resulting  dimensionless  parameters  will  be 
given  here.  The  equation  for  stresses  in  the  solid  under  equilibrium 
conditions  is: 
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the  constitutive  law  for  linear  viscoelastic  materials  is: 
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where  the  relaxation  modulus  in  shear,  G^,  can  be  expressed  as 
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Defining  the  dimensionless  parameters  as  follows 
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The  reference  quantities  selected  are  given  below  (see  Ref.  1  for 
details) 
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After  nondimensionalization,  the  equilibrium  equation  for  stresses  in 
solid  propellant  becomes: 
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The  constitutive  law  for  linear  viscoelastic  material  becomes 
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Finally,  the  complete  set  of  dimensionless  parameters  can  be  put  into 
the  following  form 
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EXPERIMENTAL  WORK 


CRACK  TIP  IGNITION  STUDY 

In  the  last  annual  report, 1  it  was  noted  that  a  bright  spot  near  the 
tip  region  of  the  propellant  crack  was  observed  before  the  convective 
ignition  front  reached  the  tip.  The  existence  of  this  spot  can  be 
attributed  to  one  of  the  following  conditions. 

1.  Ignitiov.  cf  combustible  species  evaporated  from  the  propellant 
surface  near  the  crack  entrance  region  and  transported  to  the  propellant 
tip  by  the  initial  pressure  gradient. 

2.  Ignition  or  combustion  of  the  unreacted  species  from  the 
Igniter  system  carried  to  the  crack  tip  by  the  initial  pressure  gradient. 

3.  Ignition  of  the  propellant  tip  region  by  a  mechanism  other 
than  items  1  or  2  mentioned  above. 

In  order  to  experimentally  investigate  the  phenomena  of  crack  tip 
ignition,  an  inert  crack  was  formed,  as  shown  in  Figure  11.  It  was 
designed  to  have  inert  surfaces  along  the  crack,  with  a  propellant 
sample  only  at  the  tip.  Some  modifications  were  made  in  the  combustion 
chamber  to  accommodate  the  desired  geometry.  The  crack  was  formed 
between  a  thick  plexiglass  window  and  the  bottom  surface  of  the  crack 
combustion  chamber.  To  examine  any  material  loss  of  the  propellant  at 
tha  tip,  some  go/no-go  tests  were  conducted  for  such  a  configuration. 

In  all  these  tests  the  propellant  sample  burned  and  mass  loss  were 
observed.  Since  the  tests  were  done  for  the  inert  crack,  there  was  no 
possibility  of  any  combustible  species  evaporating  from  the  crack  entrance 
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FIGURE  11.  Schematic  Diagram  of  the  Test  Rig 
for  the  Study  of  Propellant  Ignition  at  the 
Closed  End  of  an  Inert  Crack. 
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and  then  reacting  at  the  tip.  Therefore,  it  was  possible  to  rule  out 
the  first  mechanism  mentioned  (item  1).  The  tests  were  also  conducted 
with  inert  material,  e.g.,  silicon  rubber  or  clay,  instead  of  propellant 
at  the  tip  of  the  inert  crack.  High  speed  motion  pictures  of  these 
tests  revealed  no  bright  spot  near  the  tip  region,  thereby  eliminating 
the  possibility  of  the  reaction  of  some  of  the  combustible  species  from 
the  igniter  system. 

These  observations  confirm  the  fact  that  the  previously  observed 
and  reported  bright  spot  in  a  propellant  crack  is  not  caused  by  further 
reaction  of  fuel  and  oxidizer  species  resulting  from  the  incomplete 
combustion  of  the  propellant  near  the  entrance  region  of  the  crack  or 
from  the  igniter  system.  The  bright  spot  represents  the  actual  ignition 
of  the  propellant  at  the  crack  tip.  As  explained  in  the  section  entitled 
Theoretical  Work,  the  mechanism  of'  tip  ignition  appears  to  be  caused  by 
the  enhanced  transport  of  energy  in  the  gas  phase  which  is  driven  by 
strong  compression  waves. 

Numerous  experimental  test  firings  were  conducted  to  obtain  the 
effect  of  pressurization  rates  on  the  ignition  delay  time.  Propellant  A 
was  used  as  the  test  propellant  in  all  tip-ignition  tests.  A  schematic 
diagram  of  the  combustion  chamber  is  shown  in  Figure  11.  The  rectangular 
cross-sectional  channel  simulating  the  crack  was  placed  normal  to  the 
flow  direction  in  the  chamber.  A  pair  of  brass  strips  was  placed  be¬ 
tween  the  bottom  surface  of  the  combustion  chamber  and  the  innermost 
plexiglass  window  to  form  the  crack.  The  gap  width  of  the  crack  could 
be  varied  by  using  brass  strips  of  different  thicknesses.  For  all 
results  reported  here,  the  same  pair  of  brass  legs  were  used.  The  gap 
width  of  the  crack  thus  formed  was  0.051  in.  (1.3  mm).  The  length  of 
the  crack  was  about  7.48  in.  (100  mm)  and  the  width  was  0.37  in.  (9.4 
mm).  Use  of  the  plexiglass  window  as  one  of  the  surfaces  of  the  channel 
facilitated  an  unobstructed  view  of  the  ignition  phenomena.  The  chamber 
was  completely  sealed,  except  for  the  interchangeable  exit  nozzle  through 
which  the  product  gases  were  discharged  to  the  atmosphere.  The  chamber 
was  designed  for  a  maximum  static  pressure  of  800  atm.  In  this  study, 
however,  the  maximum  pressure  was  only  about  200  atm. 

A  solid  propellant  igniter  system  was  used  as  a  source  for  hot  gas 
generation.  Percussion  primers  were  used  as  the  initiators  for  ignition 
of  a  propellant  strand.  A  solenoid-activated  firing  pin  was  used  to 
trigger  the  primer.  When  the  primer  is  triggered,  the  hot  gases  flow 
over  the  propellant  strand  and  ignite  it.  The  product  gases  flow  through 
a  multi-perforated  converging  nozzle  into  the  main  chamber.  The  pressur¬ 
ization  rate  of  the  chamber  can  be  varied  by  altering  the  dimesions  of 
the  propellant  strip  in  the  igniter,  by  changing  the  dimensions  of  the 
multiperforated  nozzle,  or  by  changing  the  exit  nozzle  <j>f  the  chamber. 

In  the  present  study,  pressurization  rates  up  to  8  x  10  atm/s  were 
obtained  at  the  crack  tip. 
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A  block  diagram  of  the  remotely  controlled  ignition  and  high-speed 
photography  system  is  shown  in  Figure  12.  The  solenoid  which  triggers 
the  percussion  primer  is  activated  by  an  event-switch  built  into  the 
high-speed  camera.  After  the  camera  is  switched  on,  a  few  milliseconds 
are  required  to  reach  the  desired  framing  rate.  This  delay  time  can  be 
related  to  the  footage  of  the  film  that  passes  through  the  camera  before 
the  attainment  of  the  prescribed  framing  rate.  The  event,  switch  of  the 
camera  is  controlled  by  a  footage  indicator.  During  the  operation,  the 
footage  indicator  is  set  to  a  preselected  value  corresponding  to  the 
desired  framing  rate.  After  a  certain  length  of  film  (preset  on  the 
footage  indicator)  runs  through  the  camera,  the  event  switch  closes, 
thereby  activating  the  relay.  The  relay,  in  turn,  activates  two  other 
switches;  one  triggers  the  solenoid  to  initiate  the  ignition  event,  and 
the  other  closes  the  common-time  switch  on  a  light  emitting  diode  (LED) 
driver  unit.  At  the  instant  when  the  common-time  switch  closes,  the  LED 
driver  unit  generates  a  2  ms  pulse.  The  pulse  is  recorded  on  a  magnetic 
tape  and  is  simultaneously  marked  on  the  film.  These  common-time  marks 
are  used  for  the  time  correlation  of  the  data  recorded  on  the  tape  and 
the  film. 

The  data  acquisition  system  consisted  of  two  major  parts:  (1)  the 
pressure  recording  system,  and  (2)  the  ignition  event  recording  system. 

A  block  diagram  of  the  data  acquisition  system  is  shown  in  Figure  13. 
Pressure  measurements  were  taken  at  four  locations  along  the  crack,  in¬ 
cluding  one  at  the  crack  entrance  and  another  at  the  tip.  Piezoelectric 
quartz  transducers  (Kistler  model  601  B)  were  used  to  measure  the  pres¬ 
sure.  These  transducers  have  a  rise  time  of  1.5  us,  natural  frequency 
of  300  kHz,  and  can  accurately  record  pressures  up  to  18,000  psi  (1225 
atm).  The  transducers  were  mounted  in  a  water-cocled  adapter  (Kistler 
model  628C)  which  prevented  drifting  or  damage  caused  by  excessive  heat. 
The  signals  from  the  pressure  transducers  were  carried  through  an  in¬ 
sulated,  high-impedance,  coaxial  cable  [Model  121  M(x) ]  with  a  capaci¬ 
tance  of  30  pF/foot,  to  a  charge  amplifier  (Kistler  model  504E) .  The 
amplified  signals  were  recorded  on  FM  channels  of  a  14-channel  AM/FM 
tape  recorder  (Hewlett-Packard  model  3924E). 

Two  separate  systems  were  used  to  simultaneously  record  the  igni¬ 
tion  event.  A  high-speed,  16-mm,  motion  picture  camera  (Hycam  Model 
K20S4E-115)  was  used  to  observe  the  complete  ignition  event  in  the 
chamber.  The  camera  was  capable  of  filming  at  a  maximum  of  16,000 
pictures  per  second  and  was  equipped  with  a  dual  LED  timing  system  which 
recorded  the  common-time  and  timing  marks  on  the  film.  The  LED  driver 
unit,  a  common-time  pulse  and  high  frequency  (up  to  10  kHz)  timing 
generator,  was  used  to  operate  the  LED  timing  system.  A  photodiode 
system  was  also  used  to  observe  ignition  at  the  tip.  A  camera  lens 
(Schneider-Kreuznauch,  150  mm)  was  used  to  focus  the  tip  region  onto  a 
photodiode  (E.G.  and  G.  model  SGD040A) .  When  the  tip  region  is  ignited 
the  photodiode  gives  out  a  signal  which  is  amplified  and  recorded  on  the 
magnetic  tape. 
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FIGURE  12.  Block  Diagram  of  Remotely-Controlled 
Ignition  and  Photography  System. 


41 


NWC  TP  6278 


COMBUSTION 

CHAMBER 


PRESSURE 

TRANSDUCERS 


CHARGE  AMPLIFIER 


PHOTO  DIODE 


AMPLIFIER 


HIGH  SPEED 
MAGNETIC  TAPE 


TRANSIENT  WAVE 
FORM  DIGITIZER 


VOICE  CHANNEL 
AMPLIFIER 

42 


NWC  TP  6278 


The  data  recorded  on  tape  were  digitized  using  a  4-channel  Bioma- 
tion  Model  1015  waveform  digitizer  with  a  1024-word  memory  per  channel, 
and  a  capability  of  maximum  real-time  resolution  of  10  us .  Output  from 
the  digitizer  can  be  displayed  on  an  oscilloscope  (Tektronix  type  535 
A),  or  plotted  on  an  x-y  plotter  (Hewlett-Packard  model  7044A)  for  a 
hard  copy. 

Correlated  experimental  data  are  shown  in  Figures  14  to  18.  Even 
though  pressures  were  measured  at  four  axial  locations,  only  the  traces 
at  the  entrance  and  the  tip  are  shown  in  the  figures.  The  first  discern 
ible  pressure  rise  for  each  pressure  gauge  occurred  consecutively  from 
the  crack  entrance  to  the  crack  tip.  The  pressurization  rates  recorded 
at  various  gauge  locations  along  the  crack  indicate  that  the  pressuriza¬ 
tion  rates  decrease  consecutively  downstream  of  the  crack  entrance. 

This  attenuation  can  be  attributed  to  the  inert  nature  of  the  crack.  A 
typical  pressure-time  trace  at  the  tip  can  be  divided  into  three  regions 

(1)  the  pressurization  region  where  the  pressure  continuously  increases; 

(2)  the  constant  pressure  region;  and  (3)  the  depressurization  region. 
The  uprising  portion  of  the  pressure-time  trace  at  the  crack  tip  loca¬ 
tion  is  quite  linear.  In  most  of  the  tests,  the  ignition  occurred 
duriug  the  uprising  portion  of  the  pressure-time  curve  at  the  crack  tip 
location. 

Ignition  is  defined  in  this  report  as  the  onset  of  emission  of 
luminous  light  from  the  propellant  surface.  This  luminous  light  can  be 
detected  by  either  the  high-speed  motion  picture  camera  or  the  photo¬ 
diode.  The  high-speed  camera  was  used  for  all  the  tests;  the  photodiode 
system  was  used  only  during  the  later  tests.  It  was  noted  that  the 
photodiode  system  alone  cannot  give  a  very  reliable  signal  (see  Figure 
18).  This  is  due  to  the  fact  that  the  photodiode  picks  up  light  re¬ 
flected  from  the  walls  of  the  stainless  steel  combustion  chamber  even 
before  the  onset  of  ignition  at  the  propellant  tip.  With  a  high-speed 
motion  picture  it  was  possible  to  record  a  more  complete  ignition  event, 
thereby  reducing  the  chance  of  error. 

The  common-time  signal  is  not  shown  in  any  of  the  results  presented 
in  Figures  14  to  18.  This  is  because  there  is  a  4  to  5  ms  time  lag 
between  the  common-time  signal  and  the  first  discernible  pressure  rise 
at  the  pressure  gauge  located  at  the  crack  entrance.  Since  the  region 
of  interest  is  the  uprising  portion  of  the  pressure  traces,  that  region 
has  been  expanded,  and,  therefore,  the  common-time  signal  cannot  be 
included  in  the  same  figure.  The  time-correlated  pressure  and  ignition 
data  indicated  that  the  time  of  first  discernible  pressure  rise  coin¬ 
cides  with  the  appearance  of  hot  igniter  gases  at  the  crack  entrance. 

The  time  delay  between  the  appearance  of  the  hot  igniter  gases  at  the 
crack  entrance  and  the  ignition  at  the  tip  was  obtained  from  the  motion 
picture  recording.  The  time  lag  between  the  first  discernible  pressure 
rises  at  the  crack  entrance  and  crack  tip  was  obtained  from  the  nres- 
sure-time  trace.  The  difference  between  these  two  time  lags  is  the 
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OBSERVED  ONSET  OF  IGNITION  FROM 


HIGH  SPEED  CAMERA 


FIGURE  15.  Measured  Pressure-Time  Traces  for 
Crack  Tip  Ignition  Experiment  (DNC  Test  No.  31). 
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OBSERVED  ONSET  OF  IGNITION  FROM 


FIGURE  16.  Measured  Pressure-Time  Tracts  for 
Crack  Tip  Ignition  Experiment  (DNC  Test  'to.  33). 
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FIGURE  17.  Measured  Pressure-Time  Traces  for 
Crack  Tip  Ignition  Experiment  (DNC  Test  No.  41). 
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OBSERVED  ONSET  OF  IGNITION  FROM 


FIGURE  18.  Measured  Pressure-Time  Traces  for 
Crack  Tip  Ignition  Experiment  (DNC  Test  No.  42). 
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ignition  delay  time  for  the  propellant  at  the  tip.  This  is  identical  to 
the  definition  given  in  the  section  on  Theoretical  Work,  i.e.,  the 
ignition  delay  time  is  the  time  lag  between  the  arrival  of  the  pressure 
wave  at  the  tip  and  the  incipient  ignition.  The  plot  of  ignition  delay 
time  and  pressurization  rate  is  given  in  Figure  2  (Theoretical  Work). 
Ignition  delay  times  are  plotted  against  the  pressurization  rate  at  the 
tip  (Gage  5);  no  effort  was  made  to  relate  them  to  pressurization  rates 
at  the  crack  entrance  (Gage  2) . 


INVESTIGATION  OF  THE  ONE-DIMENSIONALITY 
OF  THE  FLAME  FRONT 

Experimental  firings  were  conducted  to  determine  the  one-dimension¬ 
ality  of  the  flame  front.  The  crack  was  formed  between  a  propellant 
slab  which  was  slightly  regressed  below  the  brass  leg  assembly1  and  the 
inner  surface  of  the  plexiglass  sacrificial  window  which  was  mounted  on 
top  of  the  leg  assembly,  i.e.,  one  side  of  the  crack  was  an  inert, 
transparent,  plexiglass  window  and  the  other  side  was  the  propellant 
slab.  With  this  type  of  configuration,  it  was  possible  to  observe  the 
flame  propagation  process  by  looking  directly  at  the  burning  propellant 
surface.  A  schematic  diagram  of  the  test  sample  used  for  this  type  of 
test  is  shown  in  Figure  19.  The  gap  width  of  the  crack  was  varied 
between  0.018  and  0.060  in.  (0.46  and  1.52  mm).  The  length  of  the 
sample  was  about  7.2  in.  (183  mm),  and  the  width  was  about  0.7  in.  (17.7 
mm)  . 


Two  pressure  traces  were  obtained  for  this  configuration:  one  at 
the  crack  entrance  and  the  other  at  the  crack  tip.  The  pressure  traces 
for  DNC  Test  No.  24  are  shown  in  Figure  20.  The  two  pressure  traces  are 
quite  similar.  The  flame  spreading  event,  for  DNC  Test  No.  24  is  depicted 
in  Figure  21.  The  filming  rate  approximately  about  12,000  pictures  per 
second.  A  detailed  discussion  of  the  flame  spreading  event  is  given  in 
Table  2.  Some  of  the  general  observations  that  can  be  made  from  these 
tests  are:  (.1)  the  hot  gases  precede  the  ignition  front  along  the  crack 
and  reach  the  tip  much  earlier  than  the  convective  ignition  front; 

(2)  the  hot  gases  are  luminous  for  a  short  time,  with  the  luminosity 
disappearing  later;  this  can  be  attributed  to  the  quenching  of  the  hot 
gases;  (3)  the  flame  front  is  nonuniform  near  the  crack  entrance,  but 
becomes  quite  uniform  along  the  crack  and  near  the  tip.  It  can  be 
concluded  that,  except  near  the  entrance  region  of  the  crack,  the  one¬ 
dimensional  assumption  for  the  flame  propagation  is  quite  appropriate. 

The  nonuniformity  of  the  flame  front  near  the  entrance  region  can  be 
attributed  to  the  nonuniform  flow  at  the  entrance  of  the  rec.tangular- 
cross-sectional  crack  channel.  These  experiments  support  the  one¬ 
dimensional  assumption  used  in  the  theoretical  model. 
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FIGURE  20.  Expanded  Short-Time  Pressure 
Traces  of  DNC  Test  No.  24. 


EFFECT  OF  PROPELLANT  CRACK  GEOMETRY 

Tests  were  continued  for  the  purpose  of  studying  the  effect  of 
propellant  crack  geometry  on  the  ignition  and  flame-spreading  phenomena. 
Since  tests  are  still  in  progress,  the  detailed  results  and  discussion 
will  be  included  in  the  next  report.  A  millmount  (an  attachment  to  the 
index  table)  has  been  designed  and  fabricated  to  manufacture  tapered 
cracks  in  the  laboratory.  Using  the  millmount,  a  maximum  taper  angle  of 
about  6  degrees  can  be  manufactured.  At  very  high  chamber  Dressuriza- 
tion  rates,  the  present  high-speed  motion  picture  recording,  system  gives 
only  a  few  useful  pictures  of  the  flame-spreading  event.  A  quarter 
frame  head  camera,  capable  of  recording  up  to  44,000  pictures  per  second, 
has  been  acquired  through  the  Pennsylvania  State  University  in  order  to 
increase  the  number  of  pictures  and  to  achieve  better  time  resolution. 
Even  at  this  high  speed,  there  is  still  the  possibility  that  some  events, 
such  as  partial  or  transient  ignition,  etc.,  may  not  be  filmed  because 
of  the  discrete  nature  of  the  framing  process.  At  the  suggestion  of 
Dr.  R.  L.  Derr,  a  streak  attachment  to  the  Hycam  camera  was  obtained  to 
alleviate  this  problem. 
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FIGURE  21.  Motion  Pictures  of  the  Flame  Spread  in 
Phenomenon  for  DN’C  Test  No.  24. 
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FIGURE  21.  (Continued.) 
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TABLE  2.  Interpretation  of  the  Motion  Pictures 
in  Figure  21  (DNC  Test  No.  24). 


Picture  No. 

Time,  ms 

Observations  and  interpretations 

1 

0 

Hot  igniter  gas  enters  the  rectangular  duct  at 
the  bottom  of  the  main  chamber 

9 

0.91 

The  crack  entrance  is  ignited.  Hot  gases  pene¬ 
trate  into  the  crack. 

15 

1.59 

The  ignition  front  is  still  near  the  entrance. 
The  hot  gases  penetrate  30%  of  the  crack. 

20 

2.15 

The  hot  gases  penetrate  more  than  30%  of  the 
crack. 

23-25 

2.50-2.72 

Hot  gases  gradually  penetrate  the  whole  length 
of  the  crack.  The  ignition  front  still  near 
the  entrance  region  of  the  crack  and  appears 
nonuniform. 

26 

2.84 

The  luminosity  of  the  hot  gases  starts  to 
decrease. 

28 

3.06 

The  hot  gas  front  recedes  to  about  70%  of  the 
crack.  The  ignition  front  moves  up. 

40 

4.42 

Hot  ga3  front  disappears.  The  ignition  front 
penetrates  about  50%  of  the  crack  length. 
Ignition  front  appears  quite  uniform. 

50,60 

5.67,6.69 

The  ignition  front  gradually  moves  up  to  about 
66%  of  the  crack.  It  appears  quite  uniform. 

80,52 

9.0,10.32 

The  ignition  front  eventually  engulfs  the 
whole  crack  length. 

SUMMARY  OF  PROGRESS  AND  CONCLUSIONS 


The  progress  made  during  the  past  12  months  in  the  study  of  combus¬ 
tion  processes  inside  solid  propellant  cracks  can  be  summarized  as 
follows: 


1.  A  theoretical  model  has 
ignition  phenomena.  It  considers 
energy  conservation  equations  for 


been  developed  to  explain  the  tip- 
one-dimensional,  unsteady  mass  and 
the  gas-phase  near  the  crack  tip. 
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Experimentally  observed  pressure-time  trace  near  the  crack-tip  region 
has  been  used  instead  of  a  gas-phase  momentum  equation.  A  transient, 
one-dimensional  heat  conduction  equation  is  considered  for  the  solid 
phase. 


2.  A  computer  program  has  been  implemented,  based  upon  the  above 
theoretical  model.  It  uses  a  stable,  implicit  scheme  for  the  solution 
of  the  finite-difference  equations.  The  program  can  predict  ignition 
delay  time,  rate  of  heat  flux  to  the  propellant  surface,  temperature 
distribution,  etc.,  as  functions  of  pressurization  rates. 

5.  Experiments  were  conducted  in  an  inert  crack  to  eliminate  some 
of  the  possible  mechanisms  for  the  ignition  at  the  tip.  Test  firings 
were  conducted  in  the  inert  crack  to  obtain  the  effect  of  pressurization 
rate  at  the  crack-tip  on  the  ignition  d^lay  time  of  the  test  propellant 
placed  at  the  crack  tip. 

4.  Both  the  theoretical  and  experimental  results  of  the  tip 
ignition  tests  indicate  that  the  ignition  delay  time  decreases  as  the 
pressurization  rate  is  increased.  Theoretically  calculated  ignition 
delay  times  are  in  good  agreement  with  experimental  data.  Based  upon 
this  agreement,  the  ignition  at  the  crack  tip  is  considered  to  be  caused 
by  the  enhanced  turbulent  energy  transport  in  the  gas  phase,  which  is 
driven  by  strong  compression  waves. 

5.  The  crack  combustion  code  (CCC)  has  been  merged  with  a  struc¬ 
tural  analysis  program  (NFAP) .  Initial  computer  runs  of  the  combined 
code  showed  that  the  solid  mechanics  and  combustion  phenomena  are  inter¬ 
dependent  . 


6.  The  crack  geometry  changes,  not  only  because  of  burning,  but 
also  because  cf  mechanical  deformation.  The  combustion  phenomena  could 
either  enhance  or  retard,  depending  upon  the  variation  in  gcemetry  along 
the  crack.  More  realistic  computation  of  the  combustion  phenomena 
inside  the  crack  should  be  based  upon  the  combined  program  (CCC/NFAP) . 

7.  The  complete  set  of  governing  equations  for  describing  the 
transient  combustion  phenomena  inside  a  solid  propellant  crack  which 
undergoes  geometry  change  due  to  both  combustion  and  mechanical  deforma¬ 
tion,  has  been  nondimensionalized.  A  set  of  governing  dimensionless 
parameters  has  been  obtained  to  aid  future  data  correlation  and  analy¬ 
sis. 


8.  Experiments  were  conducted  to  investigate  the  one-dimension¬ 
ality  of  r  .e  ignition  front.  Results  indicate  that,  except  near  the 
entrance  region  of  the  solid  propellant  crack,  the  one-dimensional 
assumption  of  the  flame  front  propagation  is  quite  appropriate. 
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9.  Experiments  were  continued  to  investigate  the  effect  of  pro¬ 
pellant  crack  geometry  on  the  ignition  and  flame-spreading  phenomena.  A 
millmount  has  been  designed  and  fabricated  to  manufacture  tapered  cracks. 
A  quarter  frame  head  camera  and  a  streak  attachment  to  the  high-speed 
camera  have  been  obtained  to  facilitate  more  accurate  visual  recording 
of  the  flame- spreading  and  ignition  events. 
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Anpendix  A 

BASIC  DEFINITIONS  AND  OPERATION  RULE  FOR  FAVRE  AVERAGING 

Definitions 

o  -  p  +  o' 

p  ®  p  +  P' 

u  ■  +  u";  x  •  ?  +  T";  h  ■  h  +  h" 


Operation  Rules 
For  any  quantity,  a. 
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DIFFICULTIES  ENCOUNTERED  IN  OBTAINING  ANALYTICAL  SOLUTIONS 


The  flnel  equation  for  the  ges  phase  (Eq.  12)  as  derived  in  section  3  is: 


3H  _  fr x  3p  .1  .  3P  .  l_  ,.  afi  v 

It  \  )o  It  dxJ  U  at  a£  (keff  ax  > 


(2-1) 


If  we  make  the  following  additional  assumptions: 

1.  P  -  Pt  (1  +  ctn) 

2.  P  - 

3.  Outside  the  thermal  boundary  layer  near  the  propellant  surface  the 

y 

compression  process  is  iesentropic:  i.e.,  £  ■  t  T'“1  "(p  )Y 


pi  Ti 


4.  k  ■  2  [kj  +  kj 

Toa 

5.  k,  -  kt  [  Tt  ] 


Then  Eq.  (2-1)  becomes: 


It  - *  €  ft ln  (<1  +  ^  < 

.  Ti  i  (1  +  ^ 

Ik  [i  +<i  +  ctV  ]  ifi 

+  2  (1+ct'1)  3x2 


(2-2) 


Introducing 

T 


1=1 


e  =  Ti  -  (1  +  ctn)Y 


1 =1 


i.e. ,  t  -  Ti  [0  +  (1  +  ctn) 

I=i 

(1  +  ctn)  Y  ] 


il  .  x  [  —  + 
at  L±  L  at 


dt 


(2-3) 


at 

ax 


t  ii 

i  ax 


a2? 

axT 


t  .  32e 
1  IxZ 


61 


NWC  TP  6278 

Substituting  thsss  into  (2-2),  ws  get: 


ii  _  _  d  n  tuY 
3t  *  3x  dt  +  ct  ' 


1=4  , 

as  iLiliWV _ l|fe. 

2  (1  +  ctH)A/Y 


(2-4) 


This  equation  can  be  transformed  into  the  form  of  a  diffusion  equation  by 
the  following  transformation: 


I  ,t 


T  2  f0  \  [l  +  (1  +  Ctn) 


i=4  - 

n'  Y  ]  (1  +  ctn)Y  dt 


X  =  xc 


1_ 

2n 


1 

2 


aQ  (1  +  ct  ') 


1 

^  Y 


(t) 


ps 


which  results  in: 

30  328 

3x  "  3X? 

I.c.  t  -  o  t  ■  o:  9  "o 
B.c.l  x  -  o  -  X  »  o:  9-6 
B.c.2  x  — *»  ;X-^  -»!  0  -+■  o 

For  the  propellant,  if  we  define 
T  ST  -  T, 

p  p  i 

Equation  (1)  of  Section  3  becomes: 


3%  1  3Tp 

ox*-  2  at 

P 


I.C.  t  ■  0 

T  -  < 

P 

B.c.l  x  ■  o 

k?44 

r  ax 

B.c.2  x  » 

0  ■  0 

1_ 
39  2H 


(2-5) 

(2-6) 


(2-7) 


(2-8) 


(2-9) 


(2-10) 


I 
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In  ord«r  co  find  a  solution  Eqs.  (2-7)  and  (2-9)  must  be  solved  simultaneously. 
A  solution  for  Eq.  (2-7)  or  (2-8)  can  be  obtained  using  Duhamel's  integral 
theorem.  However,  when  the  heat  flux  is  matched  at  the  boundary  (B.c.l  of 
Eq.  2-10),  one  obtains  an  integral  relationship,  which  does  not  seem  to  be 
analytically  solvable. 
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Appendix  C 

NFAP/CCC  PROGRAM  INPUT 


Structural  Analysis  Input  Format 


I.  HEADING 

Card  1  (  18A4  ) 

Columns  Variable 

1  -  72  HED (18) 


r 

A 


Entry 

Master  heading  information  for 
use  in  labeling  output 


Card  1 
Columns 
1  -  5 
6 

7 

8 
9 

10 

11 

12  -  15 


1315  ) 
Variable 
NUMNP 
IDOF(l) 

IDOF(2) 
IDOF (3) 
IDOF (4) 
IDOF (5) 
ID0F(6) 
NEGL 

NEGNL 


MASTER  CONTROL  CARD 

(  15,  611,  14 


Entry 

Total  number  of  nodal  points 

Master  X-translation  code 
EQ.Os  admis sable 
EQ.l:  deleted 

Master  Y-translation  code 

Master  2-translation  code 

Master  X-rotation  code 

Master  Y-rotation  code 

Master  Z-rotation  code 

Number  of  linear 
element  groups 

EQ.Os  all  elements  are  nonlinear 

Number  of  nonlinear 
element  groups 

EQ.Os  all  elements  are  linear 


G5 


1 


Card  1 
Columns 
21  -  25 


26-30 
31  -  35 
36  -  40 

|  41  -  45 

\: 

i 

.i 

46-50 

!  . 


51  -  55 
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Continued 


Variable 

HODEX 


NSTEP 

IPRI 

ISAVE 

ISTAPE 


IEQB 


NRENUM 


NPLOT 


NT  IE 


ISTAT 


Entry 

Flag  indicating 
solution  mode 
EQ.O*  data  chack  and 
plot  initial  gaomatry 
EQ.l:  axacutlon 
EQ.2:  rastart 

Number  of  solution  time  steps 

Output  printing  interval 

Restart-save  time  step  interval 
EQ.O:  no  saving  for  restart 
GE.l:  interval  of  saving  step 

Flag  indicating 
saving  solution 
on  the  tape 
EQ.O:  no 
EQ.l:  yes 

Specified  block  size 

EQ.O:  calculated  automatically 

during  the  program  run 

Flag  indicating  node 
renumbering 
to  minimize  bandwidth 
,EQ.O:  no  node  renumbering 
EQ.l:  node  renumbering 

Flag  indicating 
plotting  option 
EQ.O:  no  plot 
EQ.l:  plot  initial  mesh 
EQ.2:  Plot  geometry  at 
selected  time  steps 

Tying  degrees  of 
freedom  option 
EQ.O:  not  active 
EQ.l:  active 

Flag  indicating  static 
or  dynamic  analysis 
EQ.O:  static  analysis 
EQ.l:  dynamic  analysis 


dkft 


i 

A 

j 


3 


I 

i 


1 

I 
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Card  2  <  415,  E10.4  ) 


Columns 

Variable 

Entry 

1  -  5 

ISREF 

Number  of  time  steps 
between  reforming  effective 
stiffness  matrix 

EQ.0:  default  set  to  1 

6-10 

NUMREP 

Number  of  allowable 
stiffness  reformations  in 
each  time  step 

11  -  15 

IEQUIT 

Number  of  time  steps 
between  equilibrium  Iterations 
EQ.0:  default  set  to  1 

16  -  20 

ITEMAX 

Maximum  number  of 
equilibrium  Iterations 

EQ.0:  default  set  to  15 

21  -  30 

RTOL 

Relative  tolerance  used 
to  measure  equilibrium 
convergence  _ 

EQ.0:  default  set  to  10~ 

Card  3  ( 

2F10.5  ) 

Columns 

Vac table 

Entrx 

1  -  10 

T START 

Time  of  solution  start 

11-20 

DT 

Time  step  increment 

Card  4  ( 

315,  2F10.2,  2 OX,  3F5.0 

) 

Columns 

Variable 

Entry 

1  -  5 

KPLOT 

Plotting  option  code  to 
pick  axis  orientation 

6-10 

llAXilD 

Maximum  number  of  nodes  in  any 
element  to  be  plotted 

11  -  15 

MAXNM 

Maximum  number  of  elements  in 
any  group  to  be  plotted 

16  -  25 

VLEii 

Length  of  the  horizontal  axis 
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Card 

4  Continued 

Columns 

Variable 

Entry 

26  - 

35 

WIEN 

Length  of  the  vertical  axis 

56  - 

65 

FAC 

Scale  factor  to  magnify 
displacements  of  deformed 
geometry 

66  - 

70 

CMAG(l) 

Scele  factor  to  magnify 

X  coordinates 
(including  displacements) 

71  - 

75 

CMAG(2) 

Scale  factor  to  magnify 

Y  coordinates 
(including  displacements) 

76  - 

80 

(*MAG(3) 

Scale  factor  to  magnify 

Z  coordinates 
(Including  displacements) 


Card  5  (  1515  ) 

Columns  Variable  Entry 

1-5  NPT(l)  First  selected  solution 

step  to  be  plotted 

6-10  NPT(2)  Second  selected  solution 

step  to  be  plotted 


NPT(15)  15fck  selected  solution 

step  to  be  plotted 


Card  6  (  Blank  ) 

control  card  to  supress  NFAP  displacement  printing  routines 
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III.  NODAL  POINT  DATA 


Card  1  (  Al,  14,  Al.  14,  515,  3F10.0,  15  ) 

(Input  as  many  cards  as  required  for  all  nodes) 

Columns  Variable  Entry 

1  CT  Coordinate  system  codc. 

EQ. (blank):  cartesian  (X,Y,Z) 
EQ.X:  X-cylindrical 

2-5  N  Node  (joint)  number 

GE. 1  and  LE.  NUHNP 


6  PSF  Print  suppression  flag 

(ignored  unless  N.EQ.l) 

EQ. (blank):  no  suppression 
EQ.A:  suppress  list  of  nodal  coord. 
EQ.B:  suppress  list  of  eqtn.  numbers 
EQ.C:  both  A  and  B  above 


7-10  ID ( I , N ) 


11 

-  15 

ID(2,N) 

16 

-  20 

ID(3,N) 

21 

-  25 

ID(4,N> 

26 

-  30 

ID(5,N) 

31 

-  35 

ID(6,N) 

36 

-  45 

X(N) 

46 

-  5" 

Y(N) 

56 

-  65 

Z(N) 

66 

-  70 

KN 

X-translation  boundary  code 
EQ.O:  free 
EQ.  1:  fixed 

Y-translation  boundary  code 
Z-translation  boundary  code 
X-rotation  boundary  code 
Y-rotation  boundary  code 
Z-rotation  boundary  code 
X  coordinate 
Y  coordinate 
Z  coordinate 

Node  number  increment  for 
nodal  data  generation 
EQ.O;  no  generation 
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IV.  LOAD  CONTROL  CARDS 


Card  1  (  815  ) 

Columns 
1-5 

6-10 

11  -  15 


16  -  20 


21  -  25 


26  -  30 


31-35 
36  -  40 


Variable 

NLOAD 

NLCUR 

NP7M 

IMPACT 

NPRES 

NLCP 

NDLS 

IPGD 


Entry 


Number  of  cards  used  to 
prescribe  nodal  forces 

Total  number  of  load  curves 
(time  functions) 

Maximum  number  of  points 
used  to  describe  one  of 
the  load  curves 
GE.  NLCUR 

Impact  loading 

EQ.O:  no  impact  loading 

EQ.l:  Include  impact  loading 

Maximum  number  of  elements 
in  any  element  group 
subjected  to  pressure 
(no.  of  pressure  subgroups) 

Number  of  load  curves 
used  for  pressure  load 
EQ. NLCUR 

Number  of  pressure  load  tables 

Pressure  load  with 
geometry  dependent  flag 
EQ.O:  no 
EQ.l:  yes 


V.  INITIAL  CONDITIONS 


Card  1  (  Blank  ) 

Insert  blank  to  specify  zero  intlal  conditions 
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VI.  TWO-DIMENSIONAL  CONTINUUM  ELEMENTS 


Card 


(  3X,  Al, 


214,  4X,  14,  4X,  14,  8X,  14,  8X,  14,  4X,  414  ) 


Colugog 

4 

5-8 


9-12 


17  -  20 


25  -  28 


37  -  40 


49  -  52 


57  -  60 


61  -  64 


65  -  68 


69  -  72 


Variable 

Entry 

NPAP(l) 

Enter  "2"  for  2-D  elements 

NPAR(2) 

Number  of  2— D  elements 
in  this  group 

GE.  1 

NPAR(3) 

Type  of  nonlinear  analysis 
EQ.O  or  1:  material  nonlinear 
analysis  only 

EQ.2:  total  Lagrangian 

EQ.3:  updated  Lagrangian 

NPAR(5) 

Element  type  code 

EQ.O:  axisymmetric 

EQ.l:  plane  strain 

EQ.2:  plane  stress 

NPAR(7) 

Maximum  number  of  nodes 
to  describe  one  element 

EQ.O:  default  set  to  8 

NPAR(IO) 

Numerical  integration  order 
to  be  used  in  Gauss 
quadrature  formulae 

EQ.O:  default  set  to  2 

GE.  2  and  LE.  4 

NPAR(13) 

Number  of  stress  output 
location  tables 

EQ.O:  print  stresses  at 
integration  points 

NPAR(15) 

Material  model  numbers 

EQ.l:  linear  isotropic 

EQ.9:  linear  viscoelastic 

NEAR (16) 

Number  of  different  sets 
of  material  properties 

GE.  1 

NPAR(17) 

Number  of  constants 
per  property  set 

EQ.2:  if  NPAR(15)  .EQ.  1 
EQ.4:  if (NPAR(15)  .eq.  9 

NEAR (18) 

Dimension  of  storage  array 
required  for  element  history 
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Card  2  (  15  ) 

Columns 
1  -  5 


Variable  Entry 

N  Material  property  aet  number 


Card  3A  (  2F10.0  ) 

Use  If  linear  Isotropic  material 
Columns  Variable 

1  -  10  PROP(l,N) 

11  -  20  PROP (2, N) 


Entry 

Young's  modulus 
Poisson's  ratio 


Card  3B  (  4F10.0  ) 

Use  if  linear  viscoelastic  material 


Columns 

Variable 

Entry 

1  -  10 

PROP (1  ,N) 

Bulk  modulus 

11  -  20 

PROP (2, N) 

G.o  (first  term  in  the 
relaxation  modulus  expression) 

21  -  30 

PROP (3, N) 

(G  -  G„  ),  factor  in  the 
se8ond  terra  of  the  relaxation 
modulus  expression 

31  -  40 

PROP (4, N) 

Beta,  exponential  time  constant 
in  the  relaxation  modulus  expression 

Card  4A  (  515 

) 

Card  4A  and  4B 

needed  for  each 

element 

Columns 

Variable 

Entry 

1  -  5 

M 

Element  number 

GE.  1  and  LE.  NPAR(2) 

6-10 

IEL 

Number  of  nodes  used 
to  describe  this  element 

EQ.0:  default  set  to  NPAR(7) 
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Card  4A  Continued 


Co  lumns 
11  -  15 

16  -  20 

21  -  25 


Variable 

IPS 

MTYP 

KG 


Card  4B  (  815  ) 

Columns  Variable 

1  -  5  NOD  ( 1 ) 

6-10  NOD ( 2 ) 


36  -  40 


NOD (8) 


Entry 

Number  of  the  stress  table 

to  be  used  in  stress  calculations 

EQ.O:  no  stress  output  for  the  element 

Material  property  set  number 
assigned  to  this  element 
GE.  1  and  IE.  NPAR(16) 

Node  generation  parameter 
used  to  compute  node  numbers 
for  missing  elements  (given  on 
the  first  card  of  the  sequence) 

EQ.O:  default  set  to  1 


Entry 

Global  node  number  of 
element  node  point  1 

Global  node  number  of 
element  node  point  2 


Global  node  number  of 
element  node  point  8 


VII.  APPLIED  LOADS 


Card  1  (  15  ) 

Columns 


1  -  5 


Variable  Entry 

IPCCC  Code  for  pressure  input  source 

EQ.O:  read  pressure  load 
from  input  deck 
EQ.l:  Pressure  load  will  be 
determined  through  linkage 
to  the  CCC  program 
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Card  2A  (  215  ) 

Skip  card  set  2  If  IPCCC  .EQ.  1,  otherwise  supply 
NLCUR  sets  with  2  or  more  cards  per  set 


Columns 

Variable 

Entry 

1  -  5 

NTF 

Time  function  number 

GE.  1  and  LE.  NLCUR 

6-10 

NPTS 

Number  of  points  used 
to  describe  this  time  function 
GE.  2  and  LE.  NPTM 

Card  2B  ( 

8F10.0  ) 

Columns 

Variable 

Entry 

1  -  10 

TIMV(l) 

Time  at  point  1, 

11  -  20 

RV(1) 

Function  at  point  1,  t^ 

21  -  30 

TIMV(2) 

Time  at  point  2,  t£ 

31  -  40 

RV(2) 

Function  at  point  2,  tj 

71  -  80 

RV(4) 

Function  at  point  4,  t^ 

(next  card 

if  required) 

1  -  10 

TIMV  (5) 

Time  at  point  5,  t^ 

11  -  20 

RV(5) 

Function  at  point  5,  t^ 
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Card  3A  (  315  ) 


Repeat  card  aat  3  NDLS  time* 


Columns 

Variable 

Entry 

1  -  5 

N 

Pressure  table  ID  number 

6-10 

NT  ACE 

Element  face  on  which 
the  preasure  is  acting 

il  -  15 

LTC 

Load  type  code 

EQ.O:  Default  set  to  1 
EQ.l:  normal  pressure 

EQ.2:  hydrostatic  pressure 

Card  3B 

(  F10.0  ) 

Columns 

Variable 

Entry 

1  -  10 

• 

PCD 

Fraction  of  pressure  at 
nodes  along  pressure 
loaded  surface 

EQ.  1 

Card  4 

(  215  ) 

Columns 

Variable 

Entry 

1  -  5 

NG 

Element  group  number 

6-10 

NPRSG 

Number  of  subgroups  of 

elements  In  this  element 
group 
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Card  5  (  515  ) 


lepaat  this 

card  sat  for  as  many 

as  NPRSG  cards 

Columns 

Variable 

Entry 

1  -  5 

NELPI 

First  element  number 
for  this  subgroup 

6-10 

NELPF 

Last  element  number 
for  this  subgroup 

11  -  15 

NELINC 

Element  uumber  increment 

o 

CM 

1 

I  TABLE 

Pressure  table  ID  number 

21  -  25 

HCURP 

Load  curve  ID  number 

Card  6  (  Blank  ) 

Place  evo  blank  cards  at  end  of  data  to  signal  end  of  file 
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II.  COMPUTATIONAL  CONSTANTS 
Card  1  (  5E10.4  ) 


Columns 

Variable 

Entry 

1  -  10 

XI P 

Location  of  the  inlet 
of  the  gap 

11  -  20 

XP 

Distance  downstream  at 
which  propellant  begins 

21  -  30 

XG 

Position  at  end  of 
propellant  in  the  crack 

31  -  40 

XE 

Position  at  the  end 
of  the  crack 

41  -  50 

BX 

Body  force 

Card  2  (  6E10.4  ) 

Columns 

Variable 

Entry 

1  -  10 

W 

Average  molecular  weight 
of  all  combustion  gases 

11-20 

GAMA 

Ratio  of  specific  heats 

21  -  30 

ROPR 

Density  of  solid  propellant 

31  -  40 

FKPR 

Thermal  conductivity  of 
solid  propellant 

41  -  50 

ALPHAP 

Thermal  diffusivity  of 
solid  propellant 

51  -  60 

ROUGH 

Equivalent  relative  roughness 
of  the  port  wall 
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Card  3  (  7E10.4  ) 


Columns 

Variable 

Entry 

1  -  10 

TF 

Adiabatic  flame  temperature 

11  -  20 

TPSCRI 

Temperature  at  which 
surface  ablation  occurs 

21  -  30 

A 

Pre-exponential  factor  in  the 
non-erosive  burning  rate  law 

31  -  40 

BREXP 

Pressure  exponent  in  the 
non-eiosive  burning  rate  law 

41  -  50 

EEC 

Pre-exponential  constant  in  the 
Lenoir-Robillard  erosive 
burning  relationship 

51  -  60 

EBEX 

Factor  in  exponent  of  the 
Lenoir-Robillard  erosive 
burning  relationship 

61  -  70 

TP  SIGN 

Temperature  at  which 
ignition  of  the  solid 
propellant  occurs 

Card  4  (  4E10.4  ) 

Columns 

Variable 

Entry 

1  -  10 

TGI 

Initial  gas  temperature 
in  the  crack 

11  -  20 

PCI 

Initial  pressure  in  crack 

21  -  30 

UCI 

Initial  gas  velocity 

in  the  crack  (  to  facilitate 

the  start  of  a  numerical  solution 

31  -  40 

TP  I 

Initial  propellant  temperature 
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Card  5 
Coluana 
1  -  10 


11  -  20 

21  -  30 
31  -  40 

41  -  50 

51  -  60 

61  -  70 


Card  6 
Columns 
1  -  10 


(  7110  ) 


Variable 

NITER 


NPR1NT 


NPUN 
NS  PAVE 


NBL 


NINERT 


NSLAB 


Number  of  iterative 

pradlctor  calculation 

EQ.l:  tha  solution  is  obtained  from 

quaslllnaarlxad  inhomogeneous  tarns 

EQ.2:  Predict  twice  before 

corrector  calculation 

Printing  interval,  number  of 
time  step  calculations  after 
which  output  is  printed  out 
GE.  1 

Output  punching  interval 

Number  of  time  step 
calculations  before  space¬ 
averaging  takes  place 

Number  of  time  step 
calculations  after  each  of  which 
balance  of  the  conservation  and 
governing  equation  is  performed 
and  printed  out 

Flag  to  indicate  whether 
there  is  an  inert  wall  or  net 
EQ.0:  no  inert  wall 
GT.O:  inert  wall  present 

Number  of  propellant  slabs 
forming  the  crack 
EQ.l;  there  is  one  slab  and 
one  inert  wall 

EQ.2;  there  are  2  propellant  slabs 
forming  the  crack 


3X10,  2E10.4  ) 

Variable  Entry 

NTPPDE  Flag  to  choose  method 

of  solution  to  find 
solid  propellant  temperature 
EQ.l;  P.D.E.  solution  will  be  used 
to  solve  heat  conduction  equation 
NE.l:  Goodman's  integral  solution 
will  be  used  to  solve  surface 
temperature  variation  equation 


79 


NWC  TP  6278 


Card  6  Continued 


Columns 

Variable 

Entry 

n  -  20 

1START 

Resart  option  code 

EQ.O:  no  restart 

NE.Os  program  restart 

21  -  30 

KXRITE 

Node  number  at  which  the 
in-depth  propellant  temperature 
profile  will  be  printed 
if  NTPPDE  EQ.  1* 

Print  out  occurs  after 
every  NBL  time  steps 

31  -  40 

FAC 

Variable  used  to  generate 
unequally  spaced  in— depth  grid 
EQ.O:  produces  uniform  grid 

EQ.  0.3  typically 

41  -  50 

YTPDIP 

In-depth  location  where  the 
dipstick  ignition 
criterion  applies 

Card  7  (  7E10.4  ) 

Columns 

Variable 

Entry 

1  -  10 

DELTAX 

Spacewise  step  for  the 
numerical  calculations 

11-20 

DELTA! 

Time  step  for  the 
numerical  calculations 

21  -  30 

TI 

Time  of  solution  start 

31  -  40 

TMAX 

Maximum  time  of  a  computer 
run  or  a  restart  run 

41  -  50 

EPS LON 

Small  value  to  bypass  the 
singularity  in  the  propellant 
surface  temperature 

51  -  60 

DDRG 

Friction  coefficient 
multiplication  factor 

61  -  70 

DDHC 

Heat  transfer  coeeficient 
multiplication  factor 

i 


i 


f 


J 


80 


NWC  TP  62 7B 


Card  8  (  4U0.4  ) 


Columns 

Variable 

Entry 

1-10 

TSAVI 

Time  before  which  the  space¬ 
averaging  weight  SAVMIN  Is  ured 

11-20 

TSAV2 

Time  after  which  the  apace- 
averaging  weight  SAVMAX  is  used 

21  -  30 

SAVMIN 

Space-average  percentage  weight 
of  the  values  at  the  neighboring 
mesh  points  specified  for  the  time 
before  TSAVi 

31  -  40 

SAVMAX 

Space-average  percentage  weight 
of  the  values  at  the  neighboring 
mesh  points  specified  after  TSAV2. 
Between  TSAVI  and  RSAV2  there  is  a 
linear  increase  in  the 
space-average  weight 

Card  9  (  4E10.4  ) 

Column? 

Variable 

Entry 

1  -  10 

TTT1 

Time  before  which  the 
Crank-Nicolson  parameter  has 
the  value  THETAN 

11  -  20 

TTT2 

Time  after  which  the 

Crank-Nicolson  parameter  has 
the  value  THETAX 

21  -  30 

THETAN 

Value  of  the  Crank-Ntcolsun 
parameter  Tor  time  before  TTT1 

31  -  40 

THETAX 

V"iue  of  tne  Crank-Nicolson 
parameter  for  time  after  TTT2. 
Between  TTT1  and  TTT2,  a  linear 
increase  is  specified 
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III.  GEOMETRY  OP  CRACK  SAMPLE 


Card  1  (  15  ) 


Columns 

Variable 

1  -  5 

NAPDVX 

Card  2  (  3F15.6  ) 

Supply  NAPDVX  sets 

of  card  2 

Columns 

Variable 

1  -  15 

XAP 

16  -  30 

BPD 

31  -  45 

APD 

Card  3  <  3F15.6  ) 

Supply  NAPDVX  secs 

of  card  3 

Columns 

Variable 

1  -  15 

XAP 

16  -  30 

VFD 

31  -  45 

THETAG 

Entry 

Number  of  points  in  the 
Input  table  of  crack  geometry 


Entry 

Distance  from  the 
entrance  of  the  crack 

Burning  perimeter  of 
the  propellant  crack 

Port  area  of  the  crack 


Entry 

Distance  from  the 
entrance  of  the  crack 

Total  wetted  perimeter 
of  the  propellant  crack 

Angle  between  propellant 
wall  and  the  longitudinal  axis 
of  the  crack 
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IV.  TEST  CONDITIONS 


Card  1  (  15  ) 

Columns 
1-5 


Variable  Entry 

NTSTP  Number  of  points  used  to 

specify  the  chamber  pressure 
end  temperature  of  rocket 
with  respect  to  time 


Card  2  (  3EI2.4  ) 

Supply  NTSTP  sets  of  card  2 


Columns 

Variable 

Entry 

1  -  12 

PC 

Chamber 

pressure  of 

rocket 

13  -  24 

TC 

Chamber 

temperature 

of  rocket 

25  -  36 

TST? 

Time 
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NOMENCLATURE 


A 

A1 

A 

CP 

P 

e 

g 

G1 

h 


c 

L 

k 

"eff 

K 

M 

w 

P 


•;  b 

'  w 
R 
Re 
T 
T„ 


fpi 

ps 

t 

u 

X 

z 

a 

B 

Y 

6 

c 

\ 

u 


Axial  transformation  constant  for  the  gaseous  region,  cm 
Axial  transformation  constant  for  the  solid  region,  cm”^ 
Cross-sectional  area  of  the  crack,  cm^ 

Specific  heat  at  constant  pressure,  cal/g-K 


-1 


g-cm/gf-s2 
cal/ cm“-s-K 

Length  of  the  crack,  cm 
Thermal  conductivity  cal/cm-s-K 
k  +  k 

C  2 
Bulk  modulus,  gf/cm 

Molecular  weight,  g  /g-mole 

.  2  1 

Pressure,  gf/cm 
Prandtl  number  ^ 

Pressure,  gf/cm 

Burning  perimeter,  cm 

Wetted  perimeter  of  the  port,  cm 


Strain  deviator 

Acceleration  of  gr  .vity,  conversion  factor. 
Relaxation  Modulus,  g^/cm^ 

Enthalpy  of  gas,  cal/g 

Local  convective  heat  transfer  coefficient. 


Specific  gas  constant  for  the  combustion  gases,  g^-c.m/g-K 
Reynolds  number  1 

Temperature  (without  subscript,  gas  temperature),  K 
Adiabatic  flame  temperature  of  the  solid  propellant,  K 
Initial  propellant  temperature,  K 
Propellant  surface  temperature,  K 


Time,  s 

Gas  velocity,  cm/s 

Axial  coordinate,  measured  from  propellant  surface,  -m 

Transformed  nondimensional  axial  coordinate 

Thermal  diffusivity,  cm^/s 

Time  constant,  s~l 

Ratio  of  specific  heats 

Gap  width  of  the  crack,  cm 

Strain  tensor 

Thermal  conductivity,  cal/cm-s-K 
Gas  viscosity,  g/cm-s 
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3 

p  Density  (without  subscript,  gas  density),  g/cm 

o  Stress  tensor  ^ 

<J>  Dissipation,  cal/cm  -s 

Subscripts 

c  Rocket  chamber 

cri  Critical  condition  for  surface  ablation 

fp  Flame  propagation 

i  Initial  value 

ign  Ignition  condition 

p,  pr  Propellant 

t  Turbulent 


Superscripts 

Dimensionless  parameters  or  Favre  averaged  quantity 
*  Reference  parameters 

—  Reynold’s  or  time  averaged  quantity 

'  Fluctuating  quantities  in  Reynold's  average  method 

"  Fluctuating  quantities  in  Favre' s  average  method 
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